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Abstract 

This  study  investigated  the  frequency  response  of  a  cantilever  beam  with  the 
intent  of  establishing  a  reliable  nondestructive  method  of  damage  detection.  The 
test  specimens  were  twelve  A1  2024  T3  beams,  each  identical  excepting  a  unique  ec¬ 
centrically  located  notch,  and  one  reference  notchless  beam.  The  machined  notches 
varied  in  length  and  location  to  simulate  varying  degrees  of  damage.  Laser  doppler 
vibrometry  enabled  the  data  acquisition.  The  changes  in  natural  frequencies  were 
correlated  to  notch  length  and  notch  location.  A  comparison  of  eccentric  and  cen¬ 
tered  notch  influence  on  the  natural  frequencies  was  also  discussed. 
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INTERNAL  DAMAGE  DETECTION  AND  ASSESSMENT  IN 
BEAMS  USING  EXPERIMENTAL  NATURAL  FREQUENCIES 

I.  Introduction 

The  present  study  investigated  experimentally  the  use  of  natural  fre¬ 
quency  as  a  global  diagnostic  parameter.  Notches  were  eccentrically  placed  in  a 
beam  to  simulate  a  delamination,  and  an  acoustic  horn  causes  the  beam  to  vibrate. 
The  off-centered  placement  exaggerated  the  effects  on  the  frequency  modes  in  a 
vibrating  beam.  The  vibration  was  measured  and  analyzed  using  a  laser  dopplcr  vi- 
brometer.  This  study  investigated  the  natural  frequency  response  of  an  off-centered 
notch  in  comparison  to  a  centered  notch.  The  history  of  studies  on  natural  frequency 
response  in  beams  follow  in  Chapter  I.  Chapter  II  briefly  overviews  classical  theory, 
as  well  as  the  theory  behind  the  computer  models  and  the  laser  dopplcr  vibrometer. 
Chapter  III  outlines  the  experimental  equipment  and  the  procedure  followed.  The 
experimental  results  are  presented  and  compared  with  the  computer  modeling  re¬ 
sults,  and  further  discussed  in  Chapter  IV.  The  conclusions  are  presented  in  Chapter 
V. 

1 . 1  Motivation 

The  need  for  an  advanced  structural  health  monitoring  system  can  be 
found  in  several  industries.  The  continual  rise  in  air  traffic  has  burdened  the  aerospace 
industry  with  demands  for  lower  aircraft  production  and  maintenance  costs  coupled 
with  operational  efficiency  and  environmentally  friendly  operations,  while  maintain¬ 
ing  the  level  of  safety.  Health  monitoring  systems  hold  the  possibility  of  reducing 
structural  weight  and  minimizing  structural  inspection  requirements.20  These  effects 
in  turn  reduce  operation  and  maintenance  costs  significantly. 
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The  demand  for  advancement  of  structural  health  monitoring  is  also  growing 
in  the  space  industry,  as  spacecraft  size,  complexity,  and  lifetime  continue  to  in¬ 
crease.  Damage  is  likely  to  occur  due  to  exposure  to  the  space  environment,  orbital 
maneuvering  loads,  spacecraft  docking  operations,  and  collisions  with  space  debris.5 
The  common  visual  inspection  is  impractical  and/or  impossible  in  many  space  ap¬ 
plications,  such  as  unmanned  spacecraft  and  spacecraft  in  high-altitude  orbits.  A 
reliable  method  of  in-orbit  damage  detection  of  both  location  and  extent  is  critical 
to  schedule  necessary  repairs  and/or  mission  changes. 

A  large  portion  of  the  costly  maintenance  required  for  the  space  shuttle  is 
examination  of  the  structural  integrity.  Fatigue,  foreign  object  debris  (FOD),  and 
operational  vibrations  cause  damage  that  is  not  always  visible  at  the  surface.  The 
present  methods  of  locating  damage  are  costly  and  time-consuming.  The  present 
study  investigates  a  possible  simpler,  less  expensive  method.  If  the  self-monitoring 
systems  in  satellites  could  detect  structural  damage,  they  could  possibly  repair  the 
damage  before  operational  ability  was  impaired. 

Global  monitoring  gives  an  immediate  evaluation  of  the  structural  health.  Vi¬ 
sual  inspections  require  extensive  time  and  labor.  Often  disassembly  is  required 
to  reach  the  important  load  bearing  structures.  Manually  performed  local  non¬ 
destructive  damage  evaluation  (NDE)  methods  are  time-consuming  and  inefficient. 
Vibration-based  global  monitoring  not  only  allows  immediate  evaluations,  but  can 
also  be  used  for  continuous  long-term  monitoring.  The  goal  would  be  to  extract  in¬ 
formation  about  the  structure’s  remaining  lifetime,  time  intervals  to  the  next  main¬ 
tenance/repair  and  immediate  shutdown.12  This  information  would  benefit  many 
holds,  including  civil  engineering  structures,  railway  systems,  wind  power  stations, 
and  others.  Global  methods  require  a  coarse  sensor  network  since  they  monitor  the 
whole  system,  but  they  are  less  sensitive  and  have  a  lower  spatial  resolution  than  the 
local  methods.  However,  these  problems  can  be  resolved  by  using  a  computational 
model.12 
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1.2  Background 

Initially,  the  concept  of  monitoring  vibration  to  determine  damage  pres¬ 
ence  concerned  rotor  cracking  in  machinery  and  railroad  equipment.  Since  then, 
vibration  testing  has  been  used  extensively  in  industry  to  monitor  rotating  machin¬ 
ery.5  The  earliest  results  concerning  these  applications  is  presented  by  Dimaroganas7 
(1970).  The  work  presented  by  Nagy,  Dousis,  and  Finch  (1978)  on  the  cracks  in  rail¬ 
road  wheels  noted  the  shift  in  resonance  frequencies  due  to  the  presence  of  cracks. 

Extensive  research  has  been  recorded  on  the  vibration  testing  of  beams  with 
included  cracks/notches.  Cawley  and  Ray4  conducted  an  experiment  to  compare  the 
natural  frequency  changes  produced  by  cracks  and  slots  in  order  to  study  the  accu¬ 
racy  of  experimentally  modeling  a  crack  with  a  slot.  They  tested  simply  supported 
steel  beams  containing  single  surface  slots  of  equal  depth  and  varying  width.  The 
report  concluded  that  even  when  the  slot  width  was  small  compared  to  the  beam 
dimensions,  the  slot  width  still  had  a  strong  influence  on  the  natural  frequency 
changes.  Therefore,  the  ability  to  detect  a  slot  does  not  necessarily  ensure  the  abil¬ 
ity  to  detect  a  crack  of  the  same  depth.  Along  the  length  of  the  beam,  the  slots 
were  located  in  the  middle  for  all  the  test  cases,  which  resulted  in  significant  natural 
frequency  changes  in  the  first  and  third  flexural  modes,  but  in  negligible  changes  in 
the  second  and  fourth  modes.  The  latter  two  modes  have  a  stress  and  displacement 
node  at  the  middle  of  the  bar.  In  general,  the  natural  frequency  changes  correlated 
well  with  reduction  in  the  mid-bar  second  moment  of  inertia  due  to  an  increase  in 
width  of  the  slot. 

Kam  and  Lee13  (1992)  formulated  a  procedure  to  identify  the  cracked  element 
and  to  estimate  the  crack  size  in  a  cantilevered  beam  from  modal  test  data.  The 
beam  was  discretized  into  a  set  of  elements  and  the  crack  was  assumed  to  be  located 
within  one  of  the  elements.  The  measured  vibration  frequencies  and  mode  shapes 
were  used  to  identify  the  location  of  the  damage  by  determining  the  cracked  element 
using  a  simple  reduced  stiffness  model.  The  crack  size  is  estimated  by  comparing 
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the  strain  energy  of  the  intact  beam  with  the  tested  beam.  The  strain  energy  of  the 
uncracked  structure  is  determined  by 


Wo=i[ro]T[Jfo][ro]  (1.1) 

where  [K0\  is  the  stiffness  matrix  of  the  uncracked  structure  and  [r0]  is  the  matrix 
of  the  nodal  displacements,  which  are  obtained  from 

M  =  [A'o]-1^]  (1.2) 


The  matrix  [F]  is  the  applied  forces  on  the  uncracked  structure.  The  strain  energy  is 
found  for  the  cracked  structure  using  the  same  method.  The  stiffness  matrix  of  the 
cracked  structure  is  derived  based  on  the  stiffness  matrix  of  the  uncracked  structure 
and  the  measured  vibration  frequencies  and  mode  shapes  of  the  cracked  structure. 
With  few  exceptions,  the  crack  depth  was  estimated  within  2%. 

Several  studies  have  investigated  the  correlation  of  damage  location  and  mode 
shapes.  Yuen22  (1984)  showed  for  a  cantilevered  beam,  the  first  mode  shape  sys¬ 
tematically  changes  with  respect  to  the  damage  location.  The  model  used  involved 
a  hollow  beam  with  a  square  cross-section.  The  effect  of  the  damage  on  the  inertia 
matrix  was  ignored.  Higher  modes  did  not  exhibit  the  same  trends.  This  deviation 
which  was  attributed  to  the  complexity  of  the  higher  mode  shapes,  and  the  intricacy 
in  normalizing  the  mode  shapes. 

Pandey,  Biswas,  and  Samman15  (1990)  studied  the  changes  in  the  curvature 
mode  shapes.  Curvature  at  a  point  is  given  by 


v 


M 

El 


(1.3) 


where  M  is  the  bending  moment,  E  is  the  modulus  of  elasticity,  and  /  is  the  second 
moment  of  inertia  of  the  cross-sectional  area.  The  curvature  was  computed  from 


1-4 


numerical  results  using  a  central  difference  approximation 

„  D+ 1  -  2  Vi  +  i 

^  h2 

where  h  is  the  length  of  the  elements.  This  study  found  that  the  absolute  difference 
between  the  curvature  mode  shapes  of  the  damaged  and  undamaged  beams  was 
largest  in  the  element  that  contained  the  notch.  Additionally,  the  absolute  difference 
increased  as  the  damage  size  increased. 

Doyle10  (1995)  used  the  spectral  element  method  to  estimate  crack  size  and 
location  in  an  aluminum  beam  containing  a  transverse  crack.  The  time  response  of 
the  beam  was  measured  using  four  strain  gages  placed  on  both  sides  of  the  crack  and 
on  both  sides  of  the  beam.  A  transverse  vibration  was  initiated,  which  resulted  in 
an  axial  stress  wave  that  propagated  along  the  beam  length  and  reflected  from  the 
crack  tips  and  the  end  of  the  beam.  This  report  models  the  beam  analytically  as  two 
beams  connected  by  a  pinned  joint  with  torsional  stiffness.  The  torsional  stiffness 
can  then  be  related  to  the  crack  size. 

Palacz,  Krawczuk,  and  Ostachowicz14  (2004)  also  used  a  spectral  finite  element 
model  and  flexural-shear  coupled  wave  propagation  in  delaminated,  multiplayer  com¬ 
posite  beams.  The  beam  free  end  velocity  response,  monitored  in  the  time  domain, 
demonstrated  large  increases  in  amplitude  as  the  stress  wave  crossed  the  crack  tips 
and  reflected  from  the  free  end.  The  times  at  which  the  stress  wave  crossed  the  crack 
tips  could  be  noted  and  used  to  derive  the  crack  length. 

For  further  information  on  vibration  testing,  Dimarogonas8  (1996)  and  Salawu19 
(1995)  provide  thorough  reviews. 

The  present  study  continues  the  studies  of  an  earlier  thesis  by  Capt  Aaron 
Reifsnyder,18  in  which  cantilevered  centerline- notched  2024  T3  and  2024  O  aluminum 
beams  were  used  to  determine  whether  natural  frequency  response  measurement  in 
beam  structures  is  a  usable  method  of  damage  detection.  The  length  of  the  beams, 


(1.4) 
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excluding  the  clamped  portion,  was  30.5  cm,  the  width  was  3.81  cm,  and  the  thickness 
of  these  beams  was  only  3.175  mm.  The  lengths  of  the  notches  were  4cm,  8cm, 
12cm,  and  16cm,  which  translated  to  13%,  26%,  39%,  and  52%  of  the  beam  length. 
Twelve  test  specimens,  each  containing  a  single  notch  corresponding  to  one  of  the 
notch  lengths  and  a  location  of  the  cantilevered  end,  the  free  end,  or  mid-bar,  were 
manufactured  plus  one  reference  notchless  beam.  However,  the  beam  containing  the 
8-cm  notch  at  the  clamped  end  was  manufactured  improperly.  Additional  beams 
containing  notches  of  lengths  20  cm  and  24  cm  were  desired  for  further  comparison. 
However,  the  residual  stresses  in  the  A1  2024  T3  effected  warping  within  beams 
containing  the  longer  notches.  Therefore,  the  20-cm  and  24-cm  notched  beams  were 
fabricated  from  A1  2024  O  along  with  another  set  of  beams  containing  notches  of 
lengths  4  cm,  8  cm,  12  cm,  and  16  cm  at  the  clamped  end,  for  comparison.  The  A1 
2024  O  beams  produced  natural  frequencies  approximately  3%  greater  than  the  A1 
2024  T3.  This  difference  was  attributed  to  the  residual  stresses  in  the  A1  2024  T3 
beams.  Each  beam  was  excited  by  a  piezoelectric  actuator,  and  the  dynamic  beam 
response  was  measured  using  a  laser  Doppler  vibrometer.  The  modal  frequencies 
and  shapes  were  characterized  for  each  test  case. 


Figure  1.1  Test  support  setup  and  amplifier. 
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Figure  1.2  Laser  vibrometer  setup. 


Figure  1.3  Beams  containing  centrally  located  notches  and  reference  notchless 
beam. 
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II.  Theoretical  Development 

The  following  sections  present  the  theory  considered.  The  first  section  describes  the 
fundamental  reasoning  behind  damage  detection  through  vibration  analysis.  The 
next  section  outlines  classical  beam  theory  and  the  assumptions  involved.  The  next 
two  sections  describe  the  composition  of  the  Cobb-Durham  and  Perel-Palazotto  finite 
element  codes  respectively.  The  final  section  discusses  the  signal  processing  of  the 
laser  vibrometer. 

2.1  Natural  Frequency  Changes 

Damage  results  in  structural  reduction  in  stiffness  and  increase  in  damping. 
A  crack  changes  the  dynamic  frequency  modes  in  a  beam  by  changing  the  stiffness 
matrix,  which  is  dependent  on  the  moment  of  inertia,  I.  The  reduction  in  stiffness 
corresponds  to  decreases  in  natural  frequencies  and  alterations  in  the  modes  of  vi¬ 
brations.  The  modal  characteristics,  i.e.  the  modal  shapes  and  modal  damping, 
shift.  The  changes  in  characteristic  properties  between  the  measured  state  and  the 
intact  state  can  be  used  to  categorize  the  extent  and  location  of  the  damage.  This 
method  of  damage  detection  requires  a  reliable  reference  model.12 

A  crack  in  a  structure  also  causes  changes  in  the  damping  because  a  plastic 
zone  appears  along  the  crack  edge.  However,  the  changes  in  damping  are  usually 
relatively  small  compared  to  the  total  damping  of  the  structure.8  The  methods  used 
for  modeling  frequency  response  to  crack  depth  are  generally  linear.  However,  for  a 
breathing  crack,  a  crack  that  opens  and  closes  during  vibration,  the  linear  constant 
proportionality  is  different  for  an  open  crack  and  a  closed  crack.  Thus,  the  crack 
will  behave  as  a  bilinear  spring. 


2-1 


2.2  Classical  Beam  Theory 


The  natural  frequencies  of  an  isotropic  cantilever  beam  can  be  determined 
using  the  analytical  approach  outlined  by  Meirovitch  (1986).  A  one-dimensional 
differential  beam  element  is  considered  with  a  shear  and  a  bending  moment  at  both 
ends,  and  an  applied  external  load  along  the  length  of  the  element. 


Figure  2.1  A  differential  beam  element. 

When  the  forces  are  balanced,  the  equation  of  motion  becomes 

[Q(x,t)  +  — dx\  -  Q(x,t )  +  f(x,t)dx  =  m(x)dx —  (2-1) 

where  x  is  the  distance  along  the  length  of  the  beam  between  0  and  L,  Q(x,t )  is 
the  transverse  shearing  force,  f(x,t)  is  the  applied  external  load,  m{x)  is  the  mass 
per  unit  length,  and  y(x,t)  is  the  transverse  displacement.  Likewise,  the  moment 
equation  of  motion  is 

[M  (x,  t)  +  —  ^'’^dx]  —M  (x,  t )  +  [Q(x,  t )  +  1 dx\dx  +  f(x ,  t)dx —  =  0  (2.2) 

ox  ox  2 
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where  M(x,  t )  is  the  reaction  moment  about  the  center  axis  of  the  beam.  If  the 
element  is  sufficiently  small,  the  dx2  terms  can  be  neglected  and  the  moment  and 
force  equations  can  be  combined  to  give 


d2M(x ,  t) 
dx 2 


+  f{x,t) 


m(x) 


d2y(x,t ) 
d  t'2 


(2.3) 


The  assumption  of  transverse  shear  negligibility  is  validated  by  considering  the  thick¬ 
ness  of  the  beam  and  notch  elements.  The  beam  is  determined  thin  enough  from 
a  comparison  of  the  beam  thickness  and  length.  The  thickness  of  the  test  beams 
of  the  present  study  is  2%  of  the  beam  length.  The  4-cm  notched  element  is  also 
considered,  as  the  most  extreme  notch  case.  Calculation  reveals  the  thickest  sub- 
laminate  is  10%  of  the  notch  length.  Therefore,  the  beam  elements  are  small  enough 
to  neglect  the  through-the-thickness  shear.  The  bending  moment  is  related  to  the 
transverse  displacement  through 

M(x,t)  =  BI(*)*0*>  (2.4) 


where  E  is  the  modulus  of  elasticity  and  I (x)  is  the  moment  of  inertia  of  the  cross- 
sectional  area.  For  a  rectangular  beam,  the  moment  of  inertia  is 


bh3 

~12 


(2.5) 


where  b  is  the  base  or  width  of  the  beam  cross-section  and  h  is  the  height.  Combining 
Equations  (2.3)  and  (2.4)  yields  the  fourth-order  differential  boundary-value  equation 
of  motion 


sVm^] +/(*,*) =™.(*%9(M) 


(2,6) 


dx1'  dx1  ‘  ’  dt1 

The  natural  frequencies  are  obtained  by  considering  the  beam  with  zero  applied 
force.  The  boundary  conditions  at  the  clamped  end  are  zero  displacement  and  zero 
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slope. 


y(0,t)  =  0 

dy\  -0 

7T“  i=0  —  U 

OX 


At  the  free  end  the  curvature  and  its  derivative  are  zero. 


(2.7) 

(2.8) 


TnTd2y. 

EI-^-Ax=i  =  0 

ox 2 


d  rr^Td 2y1 
w—  [EI-^-^\x=i  =  0 


(2.9) 


(2.10) 


dx  dx2 

The  transverse  displacement  can  be  represented  as  a  product  of  a  time-dependent 
bounded  harmonic  oscillation  with  frequency  ui,  and  an  x-dependent  displacement 
function  Y(x).  Using  the  method  of  separation  of  variables,  and  incorporating  the 
boundary  conditions,  Equation  (2.6)  reduces  to 


(2.11) 


For  a  uniform  beam  with  a  constant  mass,  modulus  of  elasticity,  and  moment  of 
inertia,  Equation  (2.11)  becomes 


ai 


dx 


- p4Y(x )  =  0 


(2.12) 


where 


4  ixjn 

El 


(2.13) 


The  general  solution  to  Equation  (2.12)  is 


Y(x)  =  C\  cosh(o;x)  +  C 2  sinh(o;x)  +  C3  cos(a jx)  +  U4  sin(o;x) 


(2.14) 


where  Ci  are  unknown  coefficients  that  are  found  by  solving  the  boundary  conditions. 
When  the  clamped  and  free  boundary  conditions  are  applied  at  x  =  0  and  x  =  L 
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respectively,  the  resulting  characteristic  equation  is 


cos (/3L)  cosh {(3L)  =  — 1 


(2.15) 


Numerically  solving  Equation  (2.15)  yields  the  resonant  eigenvalues  (3r.  In  the  ab¬ 
sence  of  a  forcing  function,  the  resonant  modes  are  equivalent  to  the  natural  modes. 
The  eigenvalues  for  the  first  eight  natural  modes  are 

Table  2.1:  Resonant  eigenvalues  for  an  isotropic  can¬ 
tilevered  beam. 


Mode 

Resonant 

Eigenvalues 

1 

1.875 

2 

4.694 

3 

7.855 

4 

10.966 

5 

14.137 

6 

17.279 

7 

20.420 

8 

23.562 

The  natural  frequencies  are  calculated  by 


^ n 


ffisqrt 


El 

rriL4 


(2.16) 


where  the  material  properties  are  in  English  units.  The  displacement  function  Equa¬ 
tion  (2.14)  can  be  rewritten  as 


sin  (/3rL)  —  sin 

(2.17) 


h  (AX) 


[(sin  (3rL— sinh  /3rL)  (sin  /3rx— sinh  /lra;)  +  (cos  (3rL+ cosh  j3rL )  (cos  (3rx— cost 
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The  first  eight  mode  shapes  are  shown  in  Figure  2.2. 


2.3  Cobb-Durham  Code 

The  multi-degree  of  freedom  system  is 

[M]x  +  [C\x  +  [K]x  =  [ F]f[t )  (2.18) 

“A  multi-degree  of  freedom  system  can  be  ideally  represented  by  a  series  of  coupled 
single  degree  of  freedom  systems.”17  The  measured  FRF  can  be  understood  as  the 
superposition  of  the  FRFs  of  the  single  degree  of  freedom  systems.  In  general  a 
system  will  vibrate  at  all  eigenfrequencies9 

flX  IT 

u(x,t)  =  ^[(Ai(72)ncos(o;nt)  +  (A2C2)  sin (unt)\  sin(— — )  (2.19) 

n=0 

where  A  and  C  are  constants  that  vary  with  n  and  n  is  the  eigenfrequency.  However, 
the  system  will  vibrate  at  a  single  frequency  when  the  initial  conditions  coincide  with 
an  eigenfunction. 

The  mode  shapes  of  a  finite  element  can  be  described  as  a  sum  of  mode  shape 
functions,  assuming  the  modes  method 

y(x)  =  ^2i>i{t)i'i{t)  (2.20) 

i= 1 

where  y(x)  is  the  transverse  displacement,  ipi(t)  is  the  mode  shape  function,  and  z/*(t) 
is  the  nodal  displacement  or  rotation  corresponding  to  one  of  the  element  degrees  of 
freedom.6  The  following  boundary  conditions  are  imposed 

Vh(0)  =  l,-0i(O)  =  ML)  =  ML)  =  o  (2.21) 

MO)  =  1,  MO)  =  ML)  =  ML)  =  0  (2.22) 


2-6 
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Figure  2.2  First  eight  modes  of  vibration  of  an  isotropic  cantilever  beam. 
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^(0)  =  1,^(0)=^(L)=^3(L)=0 
04(0)  =  1,  0(00)  =  -04  (L)  =  04  (L)  =  0 


(2.23) 


(2.24) 

The  Bernoulli-Euler  beam  theory  for  a  cantilevered  beam  describes  the  mode 
shape  with 

y(a;)=C1  +  C2(y)+C3(y)2  +  C4(y)3  (2.25) 

where  ci,  C2,  C3,  and  C4  are  constants  specihc  to  the  mode.  The  mode  shape  functions 
for  the  first  four  modes  are 


0i  =  l-3(y)2  +  2  (y)3  (2.26) 

02  =  *-2/(y)2  +  /(y)3  (2.27) 

03  =  3(y)2-2(y)3  (2.28) 

04  =  -/(y)2  +  /(y)3  (2.29) 

The  elemental  matrices  are 

kij  =  Jl  Eh^^dx  (2.30) 

rriij  =  I  pA'ipi'ipjdx  (2-31) 


The  elemental  matrices  of  each  beam  element  are  combined  to  form  the  stiffness  and 
mass  matrices  of  the  entire  structure.  The  following  matrix  is  formed 


0  I 

-M~XK  - M~lC 


(2.32) 


The  eigenvalues  and  eigenvectors  of  this  matrix  are  the  resonant  frequencies  of  the 
structure.  The  Cobb-Durham  code  assumes  a  small  damping  constant. 
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Figure  2.3  Finite  element  model  of  beam. 

The  finite  element  model  was  modified  to  include  a  notch  of  variable  length. 
The  beam  was  divided  into  35  elements.  The  notched  section  contained  6  elements 
for  every  4cm.  Three  transitional  elements  were  added  at  the  ends  of  the  notches. 
The  moment  of  inertia  in  the  elements  containing  the  notch  was  modified.  For  the 
undelaminated  portion  the  moment  of  inertia  is  described  by  Equation  (2.5).  In  the 
delaminated  portion,  the  moment  of  inertia  is 


(2.33) 


where  hu  is  the  height  of  the  upper  sublaminate  and  hi  is  the  height  of  the  lower 
sublaminate.  This  change  effectively  reduced  the  moment  of  inertia,  which  reduced 
the  stiffness  matrix  to  simulate  damage. 
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Figure  2.4  Diagram  of  sublaminate  dimensions. 

2.4  Perel-Palazotto  Code 

The  theory  and  formulation  of  the  Perel-Palazotto  code  is  presented  by  Victor 
Perel16  (2003).  This  code  only  considers  thin  beams,  which  allows  the  assumption 
of  negligible  shear  strains. 

Hamilton’s  principle  is  used  to  derive  the  differential  equations  of  motion  with 
boundary  conditions. 

5  [ * 2  J(t)dt  =  0  (2.34) 

Jtx 

where  J (t)  is  a  modified  Lagrangian  function  based  on  the  strain  energy  density  and 
the  kinetic  energy  density  of  the  system.  The  latter  two  terms  are  described  in  terms 
of  density  as 

f=^p(u2  +  w2)  (2.35) 
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(2.36) 


U  =  ^axx(u'  +  ^K)2) 

where  u  is  the  longitudinal  displacement,  and  w  is  the  transverse  displacement.  The 
nonlinear  term  \{w')2  accounts  for  the  longitudinal  force  resultants  in  the  sublami¬ 
nates. 

The  transverse  displacement  of  the  beam  is  assumed  to  have  the  form 


w(x,z,t )  =  W0(x,t)  +  Dg(x)H1(z)[Wi(x,t)  -  W0(x,t)] 


(2.37) 


where  Wo(x,  t )  is  the  transverse  displacement  at  the  center  axis  of  the  beam,  W\  (x,  t ) 
is  the  transverse  displacement  of  the  upper  sublaminate  in  the  delaminated  region 
a  <  x  <  /3.  The  first  term  gives  the  transverse  displacement  for  the  undelaminated 
portions  of  the  beam.  Dp(x)  is  a  unit  step-function,  which  makes  the  second  term 
zero  in  the  undelaminated  portions  of  the  beam. 


!1  a  <  x  <  P 

0  0  <  x  <  a  and  fi  <  x  <  L 


(2.38) 


The  Heaviside  function  influences  the  delaminated  parts  of  the  beam,  and  is  depen¬ 
dent  on  the  location  along  the  thickness  of  the  beam.  In  the  lower  sublaminate  region, 
the  Heaviside  function  zeroes  the  second  term  so  that  the  transverse  displacement 
of  the  lower  sublaminate  is  W0(x,t). 


H,(z) 


0  <  z  <  7 

1  7  <  z  <  | 


(2.39) 


Otherwise,  in  the  upper  sublaminate  region,  Wi(x,t)  is  used  for  the  transverse  dis¬ 
placement. 
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In  other  words,  the  beam  can  be  described  by  four  regions  of  displacement, 


where 


w(x,  z,  t) 


wi(x,t) 
w2  (x,t) 
w3(x,t) 
W^X, t ) 


In  the  undelaminated  part  near  the  clamped  end 
In  the  lower  sublaminate 
In  the  upper  sublaminate 
In  the  undelaminated  part  near  the  free  end 


(2.40) 


The  boundary  conditions  imposed  at  either  tip  of  the  crack  are  equal  transverse 
displacement  and  rotation  for  the  upper  sublaminate  displacement,  the  lower  sub- 
laminate  displacement,  and  the  corresponding  undelaminated  part. 

The  crack  faces  are  constrained  to  prevent  interpenetration  by  accounting  for 
a  force  that  takes  effect  when  the  two  crack  faces  touch.  This  constraint  also  uses 
the  Heaviside  function  to  zero  the  force  when  the  crack  is  open,  and  to  force  the 
difference  in  displacement  to  zero  when  the  crack  is  closed.  Analytically,  the  force 
constraint  is 


f(w  2,w3)  =  (w3-w2)(1  -  —  arctan  — — — )  =  0  (2.41) 

2  7T  £ 

where  £  is  some  small  number.  The  Lagrange  multiplier  of  this  contact  forcing 
function  is  included  in  the  modified  Lagrangian  function  J(t). 

Based  on  Euler-Bernoulli  beam  theory  and  assuming  no  longitudinal  displace¬ 
ment  along  the  middle  axis,  the  longitudinal  displacement  is 

u(x,z,t)  =  ~^~z  (2.42) 
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For  the  finite  element  formulation,  the  displacement  functions  Wo{x,t)  and 
Wi(x,t)  are  approximated  by  third  degree  Hermit  polynomials. 


where 


Wo  =  [Afi  N2  N3 


iV4] 


Wo(0)  ' 

W'(0) 

w0(l) 
Wo(0  , 


Wi  =  [Ni  N2  N3 


Na]{ 


W3( 0)  ' 
W{(0) 
w3(l) 
W[{1) 


A7  3x 2  2a;3 

^  =  i--p+-p 


„2 


2a; 


a;3 


+  p 

3x 2  2a;3 


7V-(  = 


/2 


/3 


ar 


ar’ 


n‘--j+p 


(2.43) 


(2.44) 


(2.45) 

(2.46) 

(2.47) 

(2.48) 


The  elemental  coordinates  x  and  l  correspond  to  the  x-coordinate  of  the  left  node 
of  the  finite  element  and  the  length  of  the  element  respectively.  The  displacements 
w  and  u  and  the  strain  are  expressed  in  terms  of  the  eight  nodal  parameters  6 
(derivation  not  shown  for  brevity). 


w  =  [Bw]  {6} 

(2.49) 

u  =  z[Bu]{6} 

(2.50) 

mm  + 

(2.51) 
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The  elemental  stiffness  and  mass  matrices  formulated  in  terms  of  the  nodal 


parameters  are 


rl  rn  A 

[k\  =  b  /  [B£}T -^z2[B£}dzdx 

J  0  J—h/2  £>n 

rl  rh/ 2  rl  rh/ 2 

[m]  =  b  J  J  [BU]T pz2[Bu\dzdx  +  b  J  J  [BW]T  p[Bw]dzdx 


(2.52) 

(2.53) 


where  p  and  S^l  are  material  properties.  Ultimately,  the  finite  element  formulation 
for  the  strain  energy  and  kinetic  energy  is 


v  =  \mTikm + 


(2.54) 


where  Unq  is  the  nonlinear  contribution  of  the  longitudinal  force  resultants,  and 


T  =  -{0}T  [m\{theta} 


(2.55) 


The  contact  forcing  constraint  formulated  in  terms  of  nodal  parameters  is 

fk  =  (©4fc-i  -  @4fc— 3)  (tt  -  -  arctan  €>4A'~1 - ®4k-3 )  =  0  (2.56) 

2  7T  £ 

where  k  is  a  global  node  number,  and  N  is  the  total  number  of  degrees  of  freedom 
of  the  system.  The  finite  element  system  of  equations  is 


[M]{6}  +  [/?]{©}  +  x{G]  +  {«"'}  =  {P}  (2.57) 


where 

[M]  is  the  system  mass  matrix 
[K]  is  the  system  stiffness  matrix 

{Qnl}  accounts  for  the  nonlinear  component  of  the  strain  density 

\{G }  is  a  product  of  the  penalty  function  method  imposed  on  the  constraints 

{P}  is  the  force  vector 
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The  piezoelectric  actuator  introduces  a  time-dependent  concentrated  boundary 
condition.  The  boundary  condition  involves  imposing  upon  the  zone  of  the  beam 
containing  the  actuator  a  time-dependent  bending  moment.  To  simplify,  the  bending 
moment  is  introduced  in  the  forcing  function,  rather  than  the  boundary  conditions. 
The  voltage  is  distributed  evenly  along  the  length  of  the  actuator.  The  force  imposed 
on  the  beam  is  calculated  using  the  virtual  work  principle  modified  as 

b  j  J  Sii  (2)  ~h  b  Jq  J—h/ 2  .sp  (z)  (^xx  do\  — xxdzdx  T 

+  b  Jo  f—h/2  Pc(9  +  w)5wdzdx  +  b  /0L  /^02  pp(g  +  w)6wdzdx  + 

+  b  Jo  J-h/2  pcu5udzdx  +  b  /0L  J^pi  P^uSudzdx  —  /0L  qSwdx  =  ($2.58) 

where  V  is  the  applied  voltage  across  the  piezoelectric  patch,  r  is  the  thickness  of 
the  piezoelectric  patch,  d31  is  a  matrix  element  characterizing  the  material  properties 
of  the  piezoelectric  patch,  pc  is  the  density  of  beam  material,  pp  is  the  density  of 
piezoelectric  patch  material,  and  q  is  the  unit  surface  electric  charge. 

The  finite  element  code  for  the  piezoelectric  patch  and  the  delamination  were 
integrated. 


2.5  Theory  Behind  Laser  Vibrometry 

A  laser  doppler  vibrometer  is  used  to  obtain  the  frequency  response.  For  the 
system  nix  +  cx  +  kx  =  fit ),  assuming  an  under-damped  system  (J  <  1),  the 
Lagrange  transfer  function  is 


Xjs)  1 

F(s )  s2  +  2(c vns  +  u2 


(2.59) 


The  frequency  response  function  is  then 


FRFfjj) 


Xjju) 

F(jou) 


1  +  2C(f;)-(f-)2 


(2.60) 
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The  magnitude  of  the  FRF  is  maximum  at  the  resonant  frequency 

ujr  =  LUnsqrtl  —  2(2  (2.61) 

On  lightly  damped  structures  the  resonant  frequency,  damping  frequency,  and  nat¬ 
ural  frequency  approach  the  same  frequency. 

If  the  vibration  is  periodic,  the  FRF  can  be  obtained  from  the  data  of  the 
response  and  excitation  signals.  The  FRF  at  the  available  discrete  frequencies  is  the 
ratio  of  the  response  component  to  the  input  component. 

F  RF(ui)  =  ItO  (2.62) 

F(juj) 

The  velocity  and  displacement  of  a  moving  object  can  be  determined  from  back- 
scattered  light.  The  displacement  of  the  object’s  surfaced  modulates  the  phase  of 
the  light  wave.3  Shifts  in  the  light  frequency  determines  the  velocity  of  the  objects 
surface.  The  measurement  beam  from  the  interferometer  in  the  scanning  head  is 
reflected  from  mirrors  to  a  specified  scan  point  on  the  test  specimen.  The  back 
scattered  laser  light  interferes  with  the  reference  beam  in  the  scanning  head.  A 
photo  detector  records  the  interference.  A  decoder  in  the  vibrometer  provides  a 
voltage  which  is  proportional  to  the  velocity  of  the  vibration  of  the  test  specimen 
parallel  to  the  measurement  beam.  The  voltage  is  digitized  and  processed  as  the 
vibrometer  signal. 

The  sensor  head  contains  a  heterodyne  interferometer.  With  the  help  of  a 

Bragg  cell,  the  interferometer  generates  a  high-frequency  carrier  signal.  The  carrier 

signal  is  converted  into  two  quadrature  signals,  /  and  Q,  in  a  mixing  process.  This 

quadrature  signal  pair  holds  the  same  information  as  the  carrier  signal,  but  it  is  easier 

to  process  electronically  because  it  lies  in  the  base  band.  Ideally,  the  quadrature 

signals  are  sinusoidal,  have  the  same  amplitude,  and  are  phase-shifted  90°  apart. 

This  relation  means  that  the  signals  can  be  represented  by  a  rotating  pointer  whose 
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angle  of  rotation  is  equal  to  the  interferometric  phase  difference,  A (p.  The  sense  of 
rotation  corresponds  to  the  direction  of  movement  of  the  test  specimen.  The  phase 
difference  is  calculated  by 


/( A0)  =  A  cos  A 0 

(2.63) 

Q(Acj))  =  A  sin  A0 

(2.64) 

A 0  =  arctan 

J(A0) 

(2.65) 

The  phase  difference  is  proportional  to  the  displacement  x  of  the  test  specimen 
according  to 


A  0 


47rAa; 

A 


(2.66) 


where  A  is  the  laser  wavelength.  Finally,  the  Doppler  frequency  shift  can  be  obtained 
as 


A  fit) 


4  u(t) 
A 


(2.67) 


where  u{t)  is  the  velocity  of  the  specimen.  The  phase  and  frequency  demodulators 
in  the  interferometer  recover  the  displacement  and  velocity  time  histories  from  the 
phase  and  frequencies  of  the  received  signal. 


The  processes  behind  the  laser  vibrometer  hardware  and  software  are  described 
in  the  Polytec  manual.1  Periodic  functions  can  be  described  in  a  series  as  a  sum  of 
trigonometric  functions.  Non-periodic  functions  can  be  described  as  an  integral  of 
trigonometric  functions.  The  time  signal  recorded  by  the  laser  vibrometer  can  be 
split  into  a  sum  of  vibrations  at  different  frequencies  through  a  Fourier  transforma¬ 
tion.  Recorded  data  contains  a  discrete  number  of  samples.  The  Polytec  software 
uses  the  Fast  Fourier  Transformation  to  generate  the  corresponding  frequency  spec¬ 
trum. 


The  Polytec  software  provides  built-in  alias  suppression.  Aliasing  is  a  result  of 
the  following.  Any  signal  above  the  Nyquist  frequency  will  appear  in  the  spectrum 
below  the  Nyquist  frequency  at  the  same  interval  from  that  frequency.  Therefore, 
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the  spectrum  will  appear  to  have  more  signal  lines.  For  correct  signal  processing, 
f  sample  ^  2/ 'signal,  where  f sampie  is  the  sampling  frequency,  and  fsignai  is  the  maxi- 
mum  signal  frequency.  Signals  with  frequencies  fsignai  >  Jsarfp  e  generate  artificial 

p 

spectrum  lines  at  frequencies  falias  <  Jaarfv  e .  This  corruption  of  the  spectrum  is 

p 

called  the  alias  effect.  To  avoid  aliasing,  signal  frequencies  higher  than  Jsa™ple  must 
be  suppressed.  To  reject  the  alias,  the  vibrometer  signal  and  reference  signal  are 
filtered  digitally,  requiring  oversampling,  and  filtered  analog  to  anti-alias  filters.  The 
digital  filter  attenuates  frequencies  above  the  bandwidth  with  a  sharp  cutoff,  and 
the  analog  filter  attenuates  frequencies  near  the  oversampling  frequency.  The  time 
signal  is  sampled  with  an  oversampling  frequency  which  is  four  times  bigger  than 
the  sampling  frequency.  In  the  upper  20%  of  the  frequency  response,  the  amplitude 
is  already  strongly  attenuated. 

sampling  frequency  =  2.56  Bandwidth  (2.68) 

The  digital  filters  have  a  frequency  response  near  the  oversampling  frequency  and 
near  multiples  of  the  oversampling  frequency. 

Frequency  domain  averaging  consists  of  computing  the  spectrum  from  each 
time  trace  by  means  of  an  FFT,  and  then  obtaining  the  average  spectrum  by  aver¬ 
aging  all  values  at  each  frequency.  Complex  averaging  will  reduce  noise  that  is  not 
phase  correlated  to  the  input  signal.  This  method  of  averaging  can  be  used  if  the 
phase  between  the  output  and  input  signal  is  stable.  Prior  to  complex  averaging 
the  software  relates  the  phase  of  the  vibrometer  signal  to  that  of  the  reference  sig¬ 
nal.  Complex  averaging  consists  of  adding  the  real  and  imaginary  parts  of  all  the 
measurement  values  and  then  dividing  the  result  by  the  number  of  measured  values. 

1  N 

(2'69> 

iV  71=1 
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(2.70) 


i  N  N 

5^Em)  +  *EM^)] 

JV  n= 1  n= 1 

The  coherence  function  is  a  measure  of  the  power  in  the  output  signal  caused 
by  the  input.  If  the  coherence  is  1,  then  all  the  output  power  is  caused  by  the  input. 
If  the  coherence  is  0  then  none  of  the  output  is  caused  by  the  input.  Coherence 
values  close  to  zero  indicate  a  poor  signal-to-noise  ratio,  whereas  coherence  close  to 
one  indicates  a  good  signal-to-noise  ratio.  PSV  first  calculates  the  FFTs  of  all  active 
channels.  Then  the  auto  power  spectra  (AP)  for  each  channel  and  the  cross  power 
spectra  (CP)  for  every  possible  combination  of  a  vibrometer  and  reference  channel 
are  calculated. 


APSS  =  S*  x  S 

(2.71) 

CPrv  =  R*  x  V 

(2.72) 

CPvr  =  V*  x  R  =  CP*rv 

(2.73) 

Where  /  is  the  frequency,  S  is  the  spectrum  at  any  channel,  V  is  the  spectrum  of  a 
vibrometer  channel  (output),  R  is  the  spectrum  of  a  reference  channel  (input)  and  * 
denotes  the  complex  conjugate.  Without  averaging,  the  FRF  is 

FRF - = jk = s  (274) 


With  averaging,  the  response  functions  HI  and  H2  and  the  coherence  function  COH 
are  calculated: 


rri  CPrv 

vr  ~  APrr 

(2.75) 

AP 

rrn  vv 

1  vr 

(2.76) 

COHvr  = 

ii  Zivr 

(2.77) 
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These  response  functions  do  not  depend  on  the  averaging  mode  because  CP 
signals  are  always  averaged  in  complex  averaging  mode,  and  AP  signals  are  real.  HI 
and  H2  can  be  used  to  estimate  the  FRF  of  the  measured  system.  The  phase  of  HI 
and  H2  is  equal.  HI  reflects  noise  in  the  reference  signal.  However,  H2  is  affected 
by  noise  in  the  vibrometer  signal.  With  a  large  number  of  averages,  HI  converges  to 
a  smaller  value  than  the  FRF  and  H2  converges  to  a  larger  value.  Since  the  output 
noise  is  usually  larger  than  the  input,  HI  is  usually  a  better  approximation  because 
it  is  less  sensitive  to  the  output. 
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III.  Experiment  Setup  and  Procedure 


The  experiment  was  conducted  in  the  Vibration  Laboratory  of  the  Air  Force 
Institute  of  Technology  at  Wright-Patterson  Air  Force  Base,  Dayton,  Ohio.  This 
chapter  describes  the  test  specimen,  support  equipment;,  source  of  excitation,  and 
data  acquisition  equipment  and  procedures. 

3. 1  Sample  preparation 

The  test  specimens  were  designed  to  be  similar  to  those  of  the  preced¬ 
ing  thesis  in  order  to  allow  direct  comparison.  The  material  selected  for  specimen 
manufacture  was  A1  2024-T3  because  it  was  most  abundant  in  the  machine  shop. 
In  the  earlier  thesis,  Reifsnyder  found  residual  stresses  in  the  beams  composed  of 
this  material  from  the  coldworking  process  involved  in  manufacturing.  However,  the 
residual  stresses  were  determined  to  have  minimal  effect  on  the  beam  dynamics  for 
the  notch  lengths  up  to  16cm.  The  material  properties  for  A1  2024-T3  are  a  modulus 
of  elasticity  of  73.084  GPa  and  a  density  of  2837.5  kg/m3. 

The  beams  contain  notches  of  lengths  4  cm,  8  cm,  12  cm,  and  16  cm.  The 
notches  are  located  at  the  clamped  end,  at  the  free  end,  and  in  the  middle  and  ex¬ 
tend  through  the  entire  width  of  the  beam.  The  length  of  the  beam  is  38.1  cm,  of 
which  7.62  cm  are  clamped  during  the  test,  which  leaves  30.48  cm  of  the  beam  that 
is  actually  observed.  The  width  is  3.81  cm.  Both  the  length  and  width  of  the  beam 
match  those  of  the  preceding  thesis.  However,  the  thickness  of  the  beams  containing 
the  eccentric  notches  is  0.635  cm,  twice  that  of  the  beams  containing  centralized 
notches.  The  notches  within  the  beam  were  manufactured  using  an  electrical  dis¬ 
charge  machining  (EDM)  cutting  wire,  which  required  the  drilling  of  a  1/6-cm  starter 
hole.  A  smaller  hole  greatly  increased  the  likelihood  of  snapping  the  drill  bit  off  in¬ 
side  the  specimen.  Since  the  notch  was  cut  off-centered,  halfway  between  the  center 
of  the  beam  and  one  edge,  the  starter  hole  mandated  a  thicker  beam.  The  beam 
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containing  the  4  cm  notch  located  in  the  middle  was  manufactured  incorrectly  and 
included  extra  clamp  holes  at  the  free  end. 


Figure  3.1  Beams  containing  eccentrically  located  notches. 
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Figure  3.2  Beam  containing  additional  holes. 


3-3 


The  EDM  process  is  performed  by  passing  a  ultra-high  current  through  a  thin 
wire  to  burn  away  material.  The  cutting  area  is  constantly  flooded  by  a  liquid  coolant 
(in  this  case  water)  to  avoid  thermal  conduction  from  the  wire  to  areas  of  material 
beyond  the  immediate  cut.  The  specimens  for  this  study  contained  notches  with 
widths  of  0.31  mm,  which  was  the  smallest  width  achievable  by  the  EDM  machine 
used.  Cawley  and  Ray  found  that  increasing  notch  width  had  profound  effects  on 
the  variations  in  natural  frequency.  Therefore,  they  concluded  ”  it  is  desirable  to  use 
the  narrowest  available  cutter. However,  a  thinner  notch  increases  the  probability 
of  notch  surface  contact  during  vibration. 

A  piezoelectric  transducer  was  affixed  to  one  of  the  notchless  beams.  With  one 
beam  for  each  combination  of  notch  length  and  location,  and  two  notchless  beams, 
the  total  number  of  beams  is  14. 


Figure  3.3  Beam  with  piezoelectric  actuator  affixed. 


3.2  Support  system 

The  beam  was  clamped  to  a  stainless  steel  Newport  Research  Corporation 
Model  45  cylindrical  test  stand.  The  test  stand  was  affixed  through  two  socket  head 
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cap  screws  to  a  Newport  Research  Corporation  Model  100  magnetic  base  to  provide 
a  rigid  connection.  The  magnetic  base  secured  the  test  stand  to  a  pneumatic  shaker 
table.  The  clamp  was  tightened  using  four  socket  head  cap  screws.  The  screws  were 
tightened  equally  so  that  the  beam  was  aligned  flush  with  the  clamping  plate.  In 
reality  the  beam  cannot  be  perfectly  clamped,  but  the  errors  should  be  negligible. 
The  end  of  the  beam  must  be  flush  with  the  edge  of  the  clamp.  Two  holes  in  the 
clamped  section  of  the  beam  correspond  to  two  pegs  in  the  clamp,  which  ensured 
the  proper  placement  of  the  beam  within  the  clamp. 


Figure  3.4  Test  support  setup. 
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Figure  3.5  Properly  clamped  beam  flush  with  clamp. 
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The  vertical  support  was  inspected  to  ensure  that  the  test  stand  did  not  in¬ 
troduce  vibrations  or  deflections  that  would  alter  the  measurements.  The  maximum 
deflection,  which  would  occur  at  the  upper  end  of  the  vertical  support  was  calculated 
by 


p(L) 


3EI 


(3.1) 


where  P  is  the  force,  Lcyi  is  the  height  of  the  test  stand,  E  is  the  modulus  of  elasticity, 
and  /  is  the  moment  of  inertia.  The  height  of  the  test  stand  is  15  inches.  The  force 
was  derived  from  the  acceleration  of  the  beam  measured  by  the  laser  vibrometer. 
Using  the  largest  acceleration  measured  and  the  mass  of  the  unclamped  beam,  the 
maximum  force  was  determined  to  be  255  kN.  Thus  the  maximum  end  displacement 
of  the  test  stand  was  1  mm. 


3. 3  Excitation 

The  beam  was  excited  by  an  Atlas  Sound  PD-60T  acoustic  horn.  The 
horn  was  also  affixed  to  a  Newport  Research  Corporation  Model  100  magnetic  base, 
which  eliminated  base  motion.  The  horn  was  placed  behind  the  beam,  at  a  5mm 
distance  from  the  beam  in  the  center  of  the  beam  width,  as  can  be  seen  in  Figure 
3.8.  An  off-centered  excitation  would  introduce  torsion  moments,  and  disguise  the 
bending  moments  in  the  FRF.  The  horn  excited  the  beam  in  the  center  23  mm  from 
the  clamped  end.  The  horn  was  pointed  directly  perpendicular  to  the  beam.  The 
excitation  signal  was  computed  by  the  Polytec  software  program  and  passed  through 
a  220  DN  Bogen  amplifier. 
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Figure  3.6  Acoustic  horn. 


Figure  3.7  220  DN  Bogen  amplifier. 
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Figure  3.8  Distance  of  acoustic  horn  from  beam. 
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A  notchless  case  using  a  piezoelectric  actuator  was  also  prepared.  A  description 
of  the  attachment  of  the  piezoelectric  actuator  is  available  in  Reifsnyder  (2004).  The 
piezoelectric  patch  was  placed  centered  on  the  beam  to  avoid  introducing  torsional 
moments,  at  the  same  location  as  in  the  experiments  conducted  by  Reifsnyder. 


Figure  3.9  ACX  Quick  Pack  Power  Amplifier  Model  EL  1224. 

3-4  Laser  Doppler  Vibrometer 

The  laser  Doppler  vibrometer  measures  the  velocity  of  the  beam.  The 
laser  vibrometer  system  consists  of  the  ’’TOP”,  ’’LEFT”,  and  ’’RIGHT”  Polytec 
PSV-400  Scanning  Head  supported  by  3051  Manfrotto  Tripods,  the  three  corre¬ 
sponding  Polytec  OF  V-5000  Vibrometer  Controller,  and  the  Polytec  PSV-E-400-3D 
Junction  Box.  The  scanning  head  must  be  aligned  with  the  help  of  its  built-in  level 
so  that  the  path  of  the  laser  is  parallel  to  the  floor.  The  laser  should  hit  the  beam 
as  straight  on  as  possible.  The  tripod  legs  are  adjusted  so  that  the  beam  is  level  in 
the  camera  view  of  the  laser.  The  scanning  head  is  positioned  at  the  proper  height 
by  adjusting  the  tripod. 
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The  vibrometer  controller  receives  a  high-frequency  Doppler  signal  from  the 
sensor  head.  The  initial  stage  is  a  high-frequency  stage  in  which  the  measurement 
signal  is  optimally  conditioned.  Then  in  the  decoder  stage  the  velocity  information 
is  recovered  from  the  signal.  The  final  stage,  the  low-frequency  stage  consists  of  an 
analog  filter  module  and  a  digital  interface. 

3.5  Hardware  Arrangement 

The  cylindrical  test  stand  was  tightly  attached  to  the  magnetic  base  using 
two  cap  screws.  The  clamp  pegs  were  inserted  into  the  two  holes  in  the  chosen  beam. 
The  four  cap  screws  were  tightened  equally  using  an  alien  wrench  to  provide  a  fixed 
boundary  condition.  The  stand  was  then  placed  so  that  the  largest  surface  of  the 
beam  was  facing  the  area  where  the  laser  vibrometer  setup  was  established,  and  the 
magnet  was  set  in  the  arrest  position. 

The  acoustic  horn  was  also  affixed  to  a  magnetic  base  using  two  cap  screws. 
The  horn  center  was  positioned  at  the  height  of  the  middle  of  the  beam  width.  The 
horn  was  connected  to  the  ontpnt  channel  of  the  amplifier.  For  the  test  case  in 
which  the  source  of  excitation  was  the  piezoelectric  actuator,  an  ACX  Quick  Pack 
Power  Amplifier  Model  EL  1224  was  used  to  amplify  the  generator  signal.  Switches 
on  the  amplifier  were  set  to  a  voltage  limit  of  200V  and  a  current  limit  of  200mA. 
The  amplifier  was  placed  on  a  chair  separated  from  the  table,  to  prevent  added 
vibrations  from  the  internal  cooling  fan.  Alligator  clips  leading  from  the  amplifier 
output  channel  connected  to  the  prongs  of  the  piezoelectric  actuator. 

A  signal  cable,  through  which  the  generated  signal  passed,  ran  from  the  signal 
output  channel  of  the  vibrometer  junction  box  to  the  input  channel  of  the  amplifier. 
The  generated  signal  was  also  connected  to  an  input  channel  of  the  junction  box, 
to  provide  the  reference  signal  to  the  software.  Signal  cables  ran  from  all  three 
vibrometer  control  outlets  to  input  channels  in  the  junction  box. 
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Figure  3.10  Acoustic  horn  position. 

Although  the  three-dimensional  scanning  was  not  utilized,  all  three  vibrometer 
controllers  and  scanning  heads  had  to  be  turned  on.  For  one-dimensional  measure¬ 
ments,  the  “top”  laser  vibrometer  was  used.  The  laser  vibrometer  was  connected  to 
the  rest  of  the  data  acquisition  system  through  the  junction  box.  The  laser  vibrom¬ 
eter  was  positioned  so  that  the  scanning  head  line  of  sight  was  perpendicular  to  the 
test  specimen.  The  height  of  the  tripod  was  adjusted  so  that  the  scanning  head  was 
level  with  the  center  of  the  test  specimen  and  the  laser  beam  was  perpendicular  to 
the  plane  of  the  test  specimen. 

3.6  Software  preparation 

The  software  Polytec  8.2  performed  the  data  acquisition  of  the  laser  vi¬ 
brometer.  The  following  steps  were  followed  to  prepare  the  software  for  acquisition. 
With  the  test  apparatus  centered  in  the  camera  view,  the  laser  beam  was  moved 
to  a  corner  of  the  view  to  avoid  interference  and  the  camera  was  then  focused  on 
the  test  specimen.  The  laser  beam  was  then  repositioned  on  the  beam  and  properly 
focused.  The  plane  of  the  beam  was  defined.  The  laser  was  moved  to  at  least  four 
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positions  on  the  beam  in  the  camera  view  and  the  positions  were  recorded.  These 
four  readings  by  the  laser  were  used  to  create  a  plane  of  reference.  The  scan  points 
were  then  defined  by  creating  a  grid  over  the  beam  in  the  camera  screen.  Proper 
alignment  could  be  verified  by  confirming  the  laser  beam  moved  to  the  appropriate 
position  at  the  center  of  the  scan  point  when  a  scan  point  was  selected  in  the  camera 
view. 

Each  beam  was  tested  using  one-dimensional  and  two-dimensional  grids.  The 
one-dimensional  grid  avoids  measuring  much  of  the  torsion,  but  the  two-dimensional 
grid  measures  the  torsion  and  helps  to  determine  when  the  torsion  is  interfering  with 
the  measurements.  The  one-dimensional  grid  consisted  of  one  line  of  grid  points  along 
the  center  of  the  beam.  It  was  important  to  center  the  grid  points  to  avoid  including 
measurements  due  to  torsional  moments.  The  grid  points  numbered  30  for  the  coars¬ 
est  grid.  The  grid  was  refined  up  to  62  points  when  considered  necessary  to  obtain 
an  accurate  mode  shape.  For  the  two-dimensional  grid,  three  grid  points  were  placed 
across  the  beam.  The  vibrometer  signal  level  was  observed  at  several  grid  points  to 
confirm  a  sufficient  signal  amplitude,  or  signal-to-noise  ratio,  for  measurement. 

The  following  factors  were  set  in  preparation  for  acquisition.  The  acquisition 
mode  was  Fast  Fourier  Transform  (FFT)  to  process  the  beam  velocity  with  respect 
to  frequency.  The  averaging  was  complex  and  the  software  was  set  to  remeasure  data 
points  of  poor  response  signal.  The  upper  limit  of  the  frequency  range  was  set  to  8 
kHz  or  10  kHz,  the  lowest  frequency  range  that  captured  the  8th  mode.  The  window 
is  ”  rectangle”  so  no  samples  were  removed.  Twenty  samples  were  attempted  at  each 
data  point  and  averaged.  The  number  of  samples  processed  by  the  FFT  had  to  be 
a  power  of  two.  The  velocity  of  the  vibrometer  was  set  at  2  mm/s/V.  The  tracking 
filter  used  was  fast.  The  generator,  which  controlled  the  output  of  the  acoustic  horn, 
was  pseudo-random  with  amplitude  of  IV.  For  the  test  involving  the  piezoelectric 
actuator,  the  voltage  used  was  2V.  Pseudo-random  is  a  periodic  signal,  which  has  no 
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leakage  and  is  best  for  fast  measurements  of  complete  spectra  since  all  frequencies 
are  excited  simultaneously. 

3. 7  Analysis 

The  modal  analysis  used  is  the  peak-amplitude  method  with  the  assump¬ 
tion  that  in  the  vicinity  of  a  resonance  the  total  response  is  dominated  by  the  con¬ 
tribution  of  the  mode  whose  natural  frequency  is  closest.  This  method  works  when 
the  structure  is  appropriately  damped.  If  the  structure  is  too  heavily  damped  the 
resonant  frequencies  will  be  influenced  by  more  than  one  mode,  which  violates  the 
assumption  required  for  this  method. 

The  peak-amplitude  method  of  analysis  was  performed  as  follows.  The  res¬ 
onance  peaks  were  picked  from  the  FRF  plot  and  the  frequency  of  each  peak  was 
taken  as  the  natural  frequency  of  the  corresponding  mode.  This  study  investigated 
only  the  bending  modes.  Therefore,  the  modal  shapes  at  each  peak  were  analyzed 
to  discover  whether  the  peak  represented  a  bending  mode  or  a  torsional  mode.  This 
method  assumed  that  the  mode  shapes  of  the  damaged  structure  were  close  to  the 
corresponding  mode  shapes  of  the  intact  structure.  The  classical  modal  shapes  that 
indicate  a  bending  mode  are  shown. 
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IV.  Results  and  Discussion 


4.1  Abberations  in  the  FRF  plots  and  Mode  Shapes 

The  acoustic  horn  was  originally  positioned  at  the  clamped  end,  25  mm 
from  the  clamp,  which  was  the  closest  to  the  clamp  the  acoustic  horn  could  be 
positioned  without  touching  the  support  equipment.  However,  when  the  mode  shapes 
were  analyzed,  some  of  the  beams  were  missing  a  mode  and  some  of  the  mode  shapes 
indicated  interference.  This  interference  was  thought  to  be  caused  by  the  acoustic 
horn.  The  absence  of  a  mode  may  be  due  to  the  position  of  the  acoustic  horn.  If 
the  excitation  coordinates  coincides  with  a  node  for  one  of  the  modes,  then  that 
mode  will  not  appear  as  a  resonance  on  the  FRF  plot.  Since  it  is  a  node,  the  only 
response  that  will  be  encountered  at  that  modal  frequency  will  be  due  to  the  off- 
resonant  contribution  of  all  the  other  modes.11  Further  testing  was  performed  with 
the  acoustic  horn  positioned  in  the  middle  of  the  beam  and  5  mm  from  the  free 
end.  None  of  the  modes  have  nodes  that  are  located  at  the  free  end  location  of  the 
acoustic  horn.  Therefore,  the  free  end  positioning  of  the  horn  did  not  result  in  the 
absence  of  any  resonant  peaks.  Also  the  grids  were  refined  from  30  grid  points  up 
to  62. 

The  FRF  plots  produced  by  the  acoustic  horn  in  the  three  positions  for  the 
same  beam  reveal  some  influence  of  the  excitation  position  on  the  amplitude.  For 
example,  the  beam  containing  an  8-crn  notch  at  the  clamped  end  contained  interfer¬ 
ence  in  the  eighth  mode,  as  can  be  seen  in  Figure  4.1.  The  FRF  plots  for  the  three 
acoustic  horn  positions  for  the  12-cm  clamped-end  notched  beam  are  displayed  in 
Figure  4.2. 

The  amplitude  of  the  fourth  mode  is  significantly  diminished  with  the  horn 

located  at  the  middle.  The  peaks  in  the  frequency  range  of  3000  through  6000  Hz 

are  also  reduced  in  amplitude  for  the  mid-bar  horn.  The  fifth  modal  peak  with  the 

horn  located  at  the  clamped  end  is  reduced  in  amplitude.  The  eighth  mode  was 
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better  captured  by  locating  the  horn  at  the  free  end.  The  eighth  modal  shape  with 
the  horn  located  at  the  free  end  is  displayed  in  Figure  4.3. 

A  piezoelectric  patch  was  then  attached  to  a  notchless  beam.  The  FRF  plot  of 
the  notchless  case  excited  by  the  piezoelectric  patch  revealed  clear  resonance  peaks. 

The  encircled  peaks  in  Figures  4.4  and  4.5  are  the  first  seven  modal  frequencies. 
The  peaks  of  the  acoustic  horn  FRF  are  less  clear.  More  bumps  appear  in  the  horn 
FRF  than  in  the  PZT  FRF.  The  horn  FRF  also  contains  noise  in  the  extremely  low 
amplitudes  and  in  the  frequency  region  of  7000  to  8000  Hz.  This  noise  is  common 
to  all  the  beam  cases  in  the  FRF.  The  modal  frequencies  that  fall  within  this  region 
are  difficult  to  pick  out.  The  deviation  between  the  modal  frequencies  obtained  from 
the  horn  and  the  piezoelectric  actuator  are  displayed  in  Table  4.1. 

Although  the  peaks  may  be  more  difficult  to  distinguish  for  the  acoustic  horn, 
the  modal  peaks  can  still  be  determined  by  using  the  appropriate  modal  shapes.  As 
the  table  shows,  the  resulting  modal  frequencies  of  the  notchless  beam  tested  with 
the  acoustic  horn  and  the  notchlcss  beam  tested  with  the  piezoelectric  actuator  differ 
by  less  than  2%,  and  differ  much  less  than  that  for  many  of  the  modes.  However, 
the  third  mode  deviates  significantly  more  than  the  other  modes.  Note  that  in  both 
FRFs  two  peaks  occur  around  the  third  mode.  The  mode  shapes  of  these  two  peaks 
for  the  acoustic  horn  FRF  are  similar.  The  corresponding  modal  frequencies  are 


Beam  with  12  cm  Notch  at  Clamped  End,  Mode  0  at  4942  Hz 


Figure  4.1  Eighth  modal  shape  of  beam  containing  12-cm  notch  at  clamped  end, 
excited  by  acoustic  horn  at  clamped  end. 
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Acoustic  Horn 
at  Clamped  End 

Piezoelectric 

Actuator 

Error 

52.92 

53.13 

-0.40% 

330.9 

328.1 

0.84% 

1053 

996.9 

5.33% 

1869 

1869 

0.02% 

3093 

3063 

0.96% 

4590 

4588 

0.04% 

6306 

6309 

-0.04% 

8450 

8297 

1.81% 

Table  4.1  Modal  frequencies  of  notchless  beam  with  different  excitation  sources. 


828.1  Hz  and  1053  Hz.  The  discussion  that  follows  will  show  some  large  variations 
of  the  modal  frequencies  of  the  test  cases  from  the  notchless  case  in  the  third  mode. 
This  enlarged  variation  may  be  due  to  the  possibility  that  the  third  mode  is  not 
captured  and  is  actually  between  these  two  peaks. 


FRF  of  Bar  with  12-cm  Notch  at  the  Clamped  End  with  Horn  at  Free  End 


FRF  of  Bar  with  12-cm  Notch  at  the  Clamped  End  with  Horn  at  Clamped  End 


Figure  4.2  FRF  plots  of  beam  with  12-cm  notch  located  at  the  clamped  end,  ex¬ 
cited  by  horn  at  three  locations. 
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Damping  causes  the  peaks  to  decrease  in  amplitude.  However,  the  piezoelectric 
actuator  is  affixed  to  the  test  specimen  and  will  add  a  small  damping  effect,  whereas 
the  horn  does  not  touch  the  test  specimen  and  therefore  should  not  cause  notable 
damping.  By  this  reasoning,  if  any  change  in  amplitude  is  seen,  the  peaks  of  the 
piezoelectric  actuator  FRF  plot  should  be  slightly  lower. 

The  smeared  peaks  are  more  likely  due  to  interference  from  the  excitation.  The 
acoustic  horn  may  increase  the  amplitudes  of  the  torsional  modes  if  it  is  positioned 
slightly  off-centered,  or  if  the  acoustic  vibration  from  the  mouth  of  the  horn  is  not 
symmetric.  The  mode  shapes  of  the  final  modal  frequencies  chosen  are  plotted  in 
Appendix  B. 

4-2  Changes  in  Natural  Frequency 

The  resonant  frequencies  of  the  thicker  notchless  beam  excited  by  the 
acoustic  horn  are  displayed  in  Table  4.2.  The  theoretical  results  are  calculated  using 
the  analytical  approach  briefly  described  in  Chapter  II. 

The  experimental  modal  frequencies  are  approximately  4%  to  6%  less  than 
the  theoretical  modal  frequencies,  which  means  the  experimental  beam  has  a  lower 
stiffness  than  the  theoretical  case,  which  is  expected.  The  exception  is  the  third 


Beam  with  12  cm  Notch  at  Clamped  End,  Mode  8  at  4975  Hz 


Figure  4.3  Eighth  modal  shape  of  beam  containing  12  cm  notch  at  clamped  end, 
excited  by  horn  at  free  end. 
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Theoretical 

Experimental 

Error 

56.31 

52.92 

6.02% 

352.9 

330.9 

6.24% 

988.2 

1053 

-6.55% 

1937 

1869 

3.47% 

3201 

3093 

3.38% 

4782 

4590 

4.01% 

6678 

6306 

5.57% 

8892 

8450 

4.97% 

Table  4.2  Comparison  of  theoretical  and  experimental  modal  frequencies  for  the 
notchless  case. 

mode,  which  appears  to  have  an  almost  7%  greater  stiffness.  This  aberration  is  most 
likely  due  to  the  problem  noted  above,  in  finding  the  third  mode  on  the  FRF  plot 
for  the  notchless  beam. 

The  first  eight  modal  frequencies  for  each  beam  containing  an  eccentric  notch 
are  shown  in  Table  4.3.  The  table  also  shows  the  percent  change  from  the  notchless 
case,  calculated  with  respect  to  the  notchless  frequency. 


%  change 


(Modal  frequency  of  notched  beam)  —  (Modal  frequency  of  notchless  beam) 
(Modal  frequency  of  notchless  beam) 


(4.1) 

For  each  notch  location  and  mode,  as  the  notch  length  increases,  the  modal  frequency 
decreases.  The  two  exceptions  are  the  fifth  modal  frequency  of  the  12-cm  notch 


x  100% 


FRF  of  Notchless  Bar  Excited  by  Piezoelectric  Actuator 


Figure  4.4  FRF  plot  of  notchless  beam  excited  by  piezoelectric  actuator. 
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located  in  the  middle  is  larger  than  the  fifth  modal  frequency  of  the  8-cm  notch 
located  in  the  middle,  and  the  fourth  modal  frequency  of  the  8-cm  notch  located  in 
the  middle  is  larger  than  the  fourth  modal  frequency  of  the  4-cm  notch  located  in 
the  middle.  The  modal  shapes  of  the  exceptions  are  examined  in  Figures  4.6  and 
4.7. 

The  interference  can  be  recognized  by  comparison  of  the  experimental  mode 
shapes  to  the  classical  mode  shapes.  Significant  interference  is  apparent  in  the  fourth 
and  fifth  modal  shapes  of  the  8-cm  mid-notched  beam.  The  interference  in  the  fourth 
mode  is  visible  by  the  additional  maximum  at  the  location  £  =  0.7,  whereas  in  the 
classical  fourth  mode  shape,  a  minimum  occurs.  In  the  fifth  mode,  the  first  sine  curve 
is  not  smooth,  and  the  first  two  nodes  are  moved  closer  together.  The  interference 
in  both  these  mode  shapes  occurs  at  the  notch  location.  This  interference  may  be  a 
result  of  the  notch  upper  sublaminate  and  lower  sublaminate  vibrating  at  different 
frequencies,  since  they  are  of  different  thicknesses,  and  thereby  instigating  contact. 

In  order  to  be  able  to  detect  a  4-cm  notch,  the  accuracy  of  the  measurement 
system  would  have  to  be  less  than  5%,  if  only  the  first  eight  modes  are  observed. 
Damage  is  easily  detectable  in  the  8-cm  clamped-end  notched  beam  and  in  the  12-cm 
and  16-crn  notched  beams,  where  the  natural  frequency  deviates  from  the  reference 
case  by  more  than  20%. 


FRF  of  Notchless  Bar  Exerted  by  Horn  at  the  Clamped  End 


Figure  4.5  FRF  plot  of  notchless  beam  excited  by  acoustic  horn  at  clamped  end. 
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EXPERIMENTAL  MODAL  FREQUENCIES 

Deviation  from  notchless  beam  (%) 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

52.92 

330.9 

1053 

1869 

3093 

4590 

6306 

8450 

Free  End 

4-cm 

52.81 

0.21% 

328.1 

0.84% 

1009 

4.18% 

1828 

2.21% 

2972 

3.90% 

4322 

5.84% 

6178 

2.03% 

8134 

3.74% 

8-cm 

52.50 

0.79% 

325.0 

1.77% 

987.5 

6.22% 

1731 

7.40% 

2903 

6.14% 

3966 

13.59% 

5481 

13.08% 

6881 

18.57% 

12-cm 

52.50 

0.79% 

306.3 

7.42% 

942.5 

10.49% 

1666 

10.88% 

2291 

25.92% 

2550 

44.44% 

3700 

41.33% 

5023 

40.56% 

16-cm 

47.50 

10.24% 

295.0 

10.84% 

352.5 

66.52% 

857.5 

54.13% 

1485 

51.98% 

2400 

47.71% 

2805 

55.52% 

3738 

55.76% 

Middle 

4-cm 

51.88 

1.97% 

331.3 

-0.13% 

812.5 

1.96% 

1853 

0.87% 

2887 

6.66% 

4309 

6.12% 

5844 

7.33% 

8159 

3.44% 

8-cm 

50.94 

3.74% 

328.1 

0.84% 

775.0 

6.48% 

1878 

-0.46% 

2459 

20.49% 

4013 

12.57% 

5744 

8.92% 

7703 

8.84% 

12-cm 

50.31 

4.93% 

317.5 

4.04% 

710.0 

14.33% 

1638 

12.38% 

2678 

13.41% 

3850 

16.12% 

5473 

13.21% 

7308 

13.51% 

16-cm 

50.00 

5.52% 

293.8 

11.20% 

631.3 

23.82% 

906.3 

51.52% 

985 

68.15% 

1819 

60.37% 

2865 

54.57% 

3766 

55.43% 

Clamped  End 

4-cm 

50.94 

3.74% 

318.2 

3.83% 

955.7 

9.24% 

1747 

6.60% 

2858 

7.59% 

4212 

8.24% 

6279 

0.43% 

8169 

3.33% 

8-cm 

50.63 

4.33% 

296.9 

10.27% 

871.9 

17.20% 

1747 

6.54% 

2850 

7.85% 

3760 

18.08% 

4031 

36.08% 

5725 

32.25% 

12-cm 

49.38 

6.69% 

267.5 

19.15% 

575.5 

45.35% 

880 

52.92% 

1650 

46.65% 

2593 

43.51% 

3960 

37.21% 

4975 

41.12% 

16-cm 

52.50 

0.79% 

242.5 

26.71% 

352.5 

66.52% 

857.5 

54.13% 

1485 

51.98% 

2400 

47.71% 

2805 

55.52% 

3738 

55.76% 

Table  4.3  Experimental  modal  frequencies  for  the  beams  containing  eccentrically 
located  notches. 

The  modal  frequencies  of  the  beams  containing  notches  located  at  the  free 
end  are  plotted  against  mode  number  in  Figure  4.8,  the  modal  frequencies  of  the 
beams  containing  notches  located  in  the  middle  are  plotted  in  Figure  4.9,  and  the 
modal  frequencies  of  the  beams  containing  notches  located  at  the  clamped  end  are 
plotted  in  Figure  4.10.  Polynomial  least-square  curve  fits  are  generated  and  the 
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curve-fit  equations  are  displayed  for  each  case.  A  power  least-square  curve  fit  was 
deemed  appropriate  for  the  8-crn  notch  located  at  the  clamped  since  it  contained  the 
erroneous  frequencies  and  was  not  represented  well  by  a  polynomial  equation. 


Modal  Shape  of  the  Fourth  Mode  of  the  Beam  Containing  8-cm 
Notch  Located  in  the  Middle 


Figure  4.6  Fourth  modal  shape  of  beam  containing  8-cm  notch  in  the  middle. 
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Figure  4.7 


Fifth  modal  shape  of  beam  containing  8-cm  notch  in  the  middle. 


Frequency  Modes  for  Beams  with  Notches  at  the  Free  End 
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Mode  mimhei 


Figure  4.8  Frequency  modes  of  beam  containing  eccentrically  located  notch  at  the 
free  end. 
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Frequency  Modes  of  Beams  with  Notches  in  Middle  Position 
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Figure  4.9  Frequency  modes  of  beam  containing  eccentrically  located  notch  in  the 
middle. 


Frequency  Modes  of  Beams  with  Notches  at  Clamped  End 
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Figure  4.10  Frequency  modes  of  beam  containing  eccentrically  located  notch  at 
the  clamped  end. 
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The  notch  size  may  be  determined  from  the  natural  frequency  deviation  from 
the  notchless  case.  The  beams  containing  notches  located  at  the  free  end  are  shown 
in  Figure  4.8.  The  12-cm  and  16-cm  notched  beams  are  individually  distinguishable 
at  modes  4  and  higher.  The  8-cm  notched  beam  begins  to  deviate  notably  from  the 
notchless  beam  and  the  4-cm  beam  at  mode  6.  The  notches  located  at  the  clamped 
end  display  characterizable  behavior.  The  4-cm  notched  beam  does  not  appear  to 
deviate  a  large  amount  in  natural  frequencies  from  the  notchless  case.  However, 
the  8-cm,  12-cm,  and  16-cm  notched  beams  display  a  linear  proportionality  between 
natural  frequency  deviation  from  the  notchless  case  and  length  of  notch  and  appear 
distinguishable  at  mode  4  and  higher  modes.  The  4-cm  notched  beam  here  also  does 
not  appear  to  be  discernable  from  the  notchless  beam.  The  beams  containing  notches 
in  the  middle  display  a  different  behavior.  Only  the  longest  notch  is  discernable  from 
the  other  cases.  A  notch  located  eccentrically  in  the  middle  of  the  beam  is  difficult 
to  size  properly. 

In  Table  4.4,  the  deviation  from  the  notchless  case  of  the  modal  frequencies 
calculated  from  the  curve  fits  of  the  experimental  data  are  displayed.  The  percent 
deviation  is  calculated  with  respect  to  the  notchless  case,  by  Equation  (4.1).  Modes 
four  and  higher  are  shown,  since  the  lower  modes  do  not  show  discernable  behavior 
as  stated  previously,  and  the  curve  fits  are  dominated  by  the  larger  frequencies. 

In  the  case  of  4-cm  notched  beams,  the  percent  deviation  is  small.  Therefore, 
the  differences  in  percent  deviation  between  the  different  locations  are  not  significant, 
although  the  free  end  notch  produces  a  lower  deviation  than  the  others  for  modes 
four  through  seven.  In  the  case  of  8-cm  notched  beams,  deviation  of  the  clamped-end 
notch  is  significantly  higher  than  the  free-end  and  middle  notches  in  modes  four  and 
five.  In  modes  seven  and  eight,  the  percent  deviation  of  the  clamped-end  and  the 
free-end  notches  is  equivalent,  whereas  the  percent  deviation  of  the  middle  notch  is 
significantly  lower. 


4-11 


CURVE-FIT  MODAL  FREQUENCIES 
Eccentrically  notched  beams 

Deviation  from  notchless  beam  (%) 

Mode 

4 

5 

6 

7 

8 

Notchless 

1887 

3080 

4563 

6336 

8399 

Free  End 

4-cm 

1819 

3.60% 

2974 

3.44% 

4413 

3.29% 

6136 

3.16% 

8143 

3.05% 

8-cm 

1804 

4.40% 

2799 

9.12% 

3988 

12.60% 

5371 

15.23% 

6948 

17.28% 

12-cm 

1445 

23.42% 

2124 

31.04% 

2913 

36.16% 

3812 

39.84% 

4821 

42.60% 

16-cm 

938.0 

50.29% 

1473 

52.18% 

2122 

53.50% 

2885 

54.47% 

3762 

55.21% 

Middle 

4-cm 

1724 

8.64% 

2862 

7.08% 

4292 

5.94% 

6014 

5.08% 

8028 

4.42% 

8-cm 

1611 

14.63% 

2697 

12.44% 

4073 

10.74% 

5739 

9.42% 

7695 

8.38% 

12-cm 

1578 

16.38% 

2613 

15.16% 

3910 

14.31% 

5469 

13.68% 

7290 

13.20% 

16-cm 

776.0 

58.88% 

1284 

58.31% 

1948 

57.31% 

2768 

56.31% 

3744 

55.42% 

Clamped  End 

4-cm 

1724 

8.64% 

2881 

6.46% 

4342 

4.84% 

6107 

3.61% 

8176 

2.66% 

8-cm 

1426 

24.43% 

2370 

23.06% 

3588 

21.36% 

5096 

19.56% 

6906 

17.77% 

12-cm 

1007 

46.63% 

1722 

44.09% 

2643 

42.08% 

3770 

40.50% 

5103 

39.24% 

16-cm 

931.0 

50.66% 

1469 

52.31% 

2121 

53.52% 

2887 

54.43% 

3767 

55.15% 

Table  4.4  Modal  frequencies  obtained  from  the  least-squares  curve-fit  of  the  ex¬ 
perimental  data  of  the  eccentrically  notched  beams. 

The  12-crn  notches  also  display  discernable  trends  in  the  location  of  the  notches. 
The  frequency  percent  deviations  of  the  free-end  and  clamped-end  notches  are  sig¬ 
nificantly  greater  than  the  middle  notch  percent  deviations  for  all  five  modes.  In  the 
fourth  and  fifth  modes,  the  free-end  12-crn  notch  differs  from  the  clamped-end  notch 
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by  30%  to  50%.  In  the  case  of  the  16-cm  notch,  the  middle  notch  percent  deviation 
is  greater  than  that  of  the  free-end  and  clamped-end  notches  by  more  than  5%  in 
the  fourth  and  fifth  modes. 

To  establish  the  size  of  the  notch,  the  sixth  modal  frequencies  are  investigated. 
If  the  deviation  of  the  modal  frequency  from  the  notchless  case  is  between  3%  and  6% 
then  the  notch  is  13%  of  the  beam  length.  If  the  deviation  of  the  modal  frequency  is 
between  10%  and  20%  then  the  notch  is  at  least  26%  of  the  beam  length.  The  notch 
could  also  be  39%  of  the  beam  length  if  it  is  located  in  the  middle.  If  the  deviation 
of  the  modal  frequency  is  between  35%  and  45%  from  the  notchless  case  then  the 
notch  is  39%  of  the  beam  length.  If  the  deviation  is  between  50  and  60%  then  the 
notch  is  over  50%  of  the  beam  length.  These  trends  are  presented  graphically  in 
Figures  4.11  through  4.14. 


Modal  Frequency  Percent  Deviation  of4-cm  Eccentrically 
Notched  Beams 
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Figure  4.11  Modal  frequency  percent  deviation  of  4-cm  eccentrically  notched 
beams. 
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Modal  Frequency  Percent  Deviation  of  8-cm  Eccentrically 

Notched  Beams 
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Figure  4.12  Modal  frequency  percent  deviation  of  8-cm  eccentrically  notched 
beams. 


Modal  Frequency  Percent  Deviation  of  12-cm  Notched  Beams 
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Figure  4.13  Modal  frequency  percent  deviation  of  12-cm  eccentrically  notched 
beams. 
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Modal  Frequency  Percent  Deviation  of  16-cm  Eccentrically 
Notched  Beams 
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Figure  4.14 


Modal  frequency  percent  deviation  of  16-cm  eccentrically  notched 
beams. 
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4-3  Comparison  with  Beams  Containing  Centered  Notches 

These  results  correspond  with  the  modal  frequency  trends  of  the  beams 
containing  the  centralized  notch  discussed  by  Reifsnyder.18  In  order  to  be  able  to 
compare  the  beams  of  different  thicknesses,  the  modal  frequencies  are  normalized 
using  the  corresponding  modal  frequencies  of  the  notchless  beam  of  the  appropriate 
thickness.  For  the  case  of  the  beam  containing  an  8  cm  notch  at  the  clamped  end,  the 
results  for  the  A1  2024  O  case  are  displayed  with  a  3%  correction  factor  to  account 
for  the  residual  stresses  within  the  A1  2024  T3  beams.  Again,  the  percent  deviation 
is  calculated  using  Equation  (4.1). 

The  ratio  for  the  notchless  case  presented  in  Table  4.5  is  the  eccentrically 
notched  beam  frequency  divided  by  the  centrally  notched  beam  frequency.  The 
theoretical  natural  frequency  for  a  simple  lump  mass  is 

(4.2) 

where  k  is  the  stiffness  and  m  is  the  mass  of  the  structure.  The  stiffness  is  a  function 
of  the  moment  of  inertia  /,  which  is  found  using  Equation  (2.5).  Since  the  thickness 
of  the  beam  containing  the  centered  notch  is  half  that  of  the  beam  containing  the 
eccentric  notch,  the  moment  of  inertia  of  the  former  is  one  eighth  that  of  the  latter. 
On  the  other  hand,  since  the  beams  are  isotropic  and  the  same  material,  the  mass  of 
the  centrally  notched  beam  is  half  that  of  the  eccentrally  notched  beam.  Therefore, 
the  ratio  of  the  natural  frequencies  of  the  notchless  beams  is 

—  =  sqrtf  =  2  (4.3) 

A  2 

The  experimental  ratios,  which  are  shown  for  the  notchless  case  in  Table  4.5,  are 
found  to  be  approximately  2,  as  well. 
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EXPERIMENTAL  MODAL  FREQUENCIES 

Normalized  with  notchless  beam 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

Eccentric 

52.92 

330.9 

1053 

1869 

3093 

4590 

6306 

8450 

Centered 

28.09 

172.4 

484 

945 

1566 

2342 

3268 

4359 

Ratio 

1.88 

1.92 

2.18 

1.98 

1.98 

1.96 

1.93 

1.94 

Free  End 

4-cm 

3.81x0.64x30.5  (eccentric) 

0.998 

0.992 

0.958 

0.978 

0.961 

0.942 

0.980 

0.963 

3.81x0.32x30.5  (centered) 

1.029 

1.016 

0.974 

0.986 

0.948 

0.945 

0.954 

0.952 

%  Difference 

-3.13% 

-2.50% 

-1.69% 

-0.85% 

1.30% 

-0.37% 

2.66% 

1.08% 

8-cm 

3.81x0.64x30.5  (eccentric) 

0.992 

0.982 

0.938 

0.926 

0.939 

0.864 

0.869 

0.814 

3.81x0.32x30.5  (centered) 

1.043 

0.968 

0.871 

0.820 

0.868 

0.760 

0.775 

0.798 

%  Difference 

-5.14% 

1.44% 

7.13% 

11.49% 

7.47% 

12.07% 

10.80% 

2.06% 

12-cm 

3.81x0.64x30.5  (eccentric) 

0.992 

0.926 

0.895 

0.891 

0.741 

0.556 

0.587 

0.594 

3.81x0.32x30.5  (centered) 

1.034 

0.895 

0.839 

0.697 

0.696 

0.747 

0.725 

0.699 

%  Difference 

-4.21% 

3.34% 

6.24% 

21.83% 

5.99% 

-34.37% 

-23.57% 

-17.63% 

16-cm 

3.81x0.64x30.5  (eccentric) 

0.898 

0.892 

0.335 

0.459 

0.480 

0.523 

0.445 

0.442 

3.81x0.32x30.5  (centered) 

1.006 

0.851 

0.695 

0.594 

0.668 

0.655 

0.693 

0.631 

%  Difference 

-12.08% 

4.53% 

-107.71% 

-29.55% 

-39.18% 

-25.35% 

-55.75% 

-42.60% 

Middle 

4-cm 

3.81x0.64x30.5  (eccentric) 

0.980 

1.001 

0.980 

0.991 

0.933 

0.939 

0.927 

0.966 

3.81x0.32x30.5  (centered) 

0.987 

0.989 

0.942 

0.952 

0.981 

0.844 

0.934 

0.868 

%  Difference 

-0.70% 

1.21% 

3.87% 

3.98% 

-5.08% 

10.05% 

-0.82% 

10.06% 

8-cm 

3.81x0.64x30.5  (eccentric) 

0.963 

0.992 

0.935 

1.005 

0.795 

0.874 

0.911 

0.912 

3.81x0.32x30.5  (centered) 

1.001 

0.980 

0.807 

0.768 

0.819 

0.789 

0.758 

0.799 

%  Difference 

-4.03% 

1.16% 

13.66% 

23.53% 

-3.06% 

9.81% 

16.78% 

12.35% 

12-cm 

3.81x0.64x30.5  (eccentric) 

0.951 

0.960 

0.857 

0.876 

0.866 

0.839 

0.868 

0.865 

3.81x0.32x30.5  (centered) 

0.978 

0.936 

0.756 

0.673 

0.690 

0.738 

0.728 

0.686 

%  Difference 

-2.87% 

2.50% 

11.72% 

23.21% 

20.29% 

11.97% 

16.16% 

20.67% 

16-cm 

3.81x0.64x30.5  (eccentric) 

0.945 

0.888 

0.762 

0.485 

0.318 

0.396 

0.454 

0.446 

3.81x0.32x30.5  (centered) 

0.955 

0.884 

0.671 

0.611 

0.654 

0.599 

0.607 

0.647 

%  Difference 

-1.05% 

0.45% 

11.90% 

-26.03% 

-105.42% 

-51.20% 

-33.53% 

-45.16% 

Clamped  End 

4-cm 

3.81x0.64x30.5  (eccentric) 

0.963 

0.962 

0.908 

0.935 

0.924 

0.918 

0.996 

0.967 

3.81x0.32x30.5  (centered) 

0.987 

0.985 

0.981 

1.022 

0.886 

0.887 

0.946 

0.905 

%  Difference 

-2.56% 

-2.46% 

-8.08% 

-9.35% 

4.12% 

3.34% 

4.94% 

6.36% 

8-cm 

3.81x0.64x30.5  (eccentric) 

0.957 

0.897 

0.828 

0.935 

0.922 

0.819 

0.639 

0.678 

3.81x0.32x30.5  (centered) 

0.958 

0.897 

0.973 

0.731 

0.789 

0.753 

0.795 

0.781 

%  Difference 

-0.09% 

0.07% 

-17.49% 

21.82% 

14.37% 

8.10% 

-24.45% 

-15.24% 

12-cm 

3.81x0.64x30.5  (eccentric) 

0.933 

0.808 

0.547 

0.471 

0.534 

0.565 

0.628 

0.589 

3.81x0.32x30.5  (centered) 

0.895 

0.850 

0.763 

0.723 

0.664 

0.707 

0.701 

0.676 

%  Difference 

4.12% 

-5.18% 

-39.52% 

-53.52% 

-24.42% 

-25.19% 

-11.59% 

-14.90% 

16-cm 

3.81x0.64x30.5  (eccentric) 

0.992 

0.733 

0.335 

0.459 

0.480 

0.523 

0.445 

0.442 

3.81x0.32x30.5  (centered) 

0.862 

0.795 

0.620 

0.597 

0.655 

0.619 

0.605 

0.632 

%  Difference 

13.09% 

-8.48% 

-85.31% 

-30.07% 

-36.48% 

-18.33% 

-36.02% 

-42.90% 

Table  4.5  Comparison  of  eccentrically  and  centrally  notched  beams  normalized  by 
the  notchless  cases. 

The  normalized  frequencies  deviate  less  than  10%  for  the  beams  containing 
notches  of  length  4  cm.  An  eccentrically  located  notch  and  a  centrally  located  notch 
are  comparable  for  a  small  crack.  The  normalized  frequencies  of  beams  containing 
8-cm  notches  are  also  close.  The  highest  deviations  between  the  eccentrically  and 
centrally  notched  beams  are  12%  for  the  6th  mode  of  the  notch  located  at  the  free 
end  and  24%  for  the  4th  mode  of  the  notch  located  in  the  middle. 

The  differences  become  significant  for  the  beams  containing  longer  notches. 
The  12-cm  notch  located  in  the  middle  produces  frequencies  that  are  still  within 
24%.  For  this  case,  the  eccentric  beam  actually  has  a  higher  stiffness  indicated  by 
the  larger  normalized  natural  frequencies.  The  beams  containing  12-cm  notches  at 
the  free  end  and  the  clamped  end  display  normalized  natural  frequencies  that  are 
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significantly  lower  in  the  eccentrically  notched  beam.  This  deviation  indicates  a 
lower  stiffness,  which  is  expected.  The  16-crn  notched  beams  display  this  behavior 
at  all  three  notch  locations.  The  largest  deviations  are  108%  for  the  3rd  mode  of  the 
free  end  notched  bar,  105%  for  the  5th  mode  of  the  mid-bar  notch,  and  85%  for  the 
3rd  mode  of  the  clamped  end  notched  bar.  This  suggest  that  for  longer  notches,  the 
eccentrically  placed  notch  has  a  greater  influence  on  the  natural  frequency.  However, 
the  influence  of  the  beam  thickness  differences  may  not  be  canceled  by  the  simple 
normalization. 


Modal  Frequencies  of  Beams  Containing  Centrally  Located 
Notches  at  the  Free  End 
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Figure  4.15 


Frequency  modes  of  beam  containing  centrally  located  notch  at  the 
free  end. 
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Modal  Frequencies  of  Beams  Containing  Notches  Centrally 
Located  in  the  Middle 
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Figure  4.16  Frequency  modes  of  beam  containing  centrally  located  notch  in  the 
middle. 


Modal  Frequencies  of  Beams  Containing  Centrally  Located 
Notches  at  the  Clamped  End 
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Figure  4.17  Frequency  modes  of  beam  containing  centrally  located  notch  at  the 
clamped  end. 
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Unlike  the  eccentrically  notched  beams,  the  beams  containing  centrally  located 
notches  display  a  linear  proportionality  between  notch  length  and  deviation  from  the 
reference  case  for  all  three  notch  locations.  However,  the  notches  become  discernablc 
at  a  higher  mode  than  for  the  eccentrically  notched  beams. 

In  Table  4.6,  the  deviation  from  the  notchless  case  of  the  modal  frequencies 
calculated  from  the  curve  fits  of  the  experimental  data  are  displayed.  Modes  four 
and  higher  are  shown,  to  compare  with  the  eccentrically  notched  beam  results. 

Observation  of  the  percent  deviation  from  the  notchless  case  yields  the  fol¬ 
lowing  comparison  of  the  centered  and  eccentric  notches.  The  centered  notch  has 
a  larger  deviation  from  the  notchless  beam  than  the  eccentric  notch  for  the  5th 
through  8th  modal  frequencies  for  the  beams  with  4-cm  and  8-crn  notches  at  each  of 
the  three  locations.  This  observation  also  applies  to  the  12-cm  middle  notched  beam. 
However,  the  beams  containing  the  12-cm  notch  at  the  free  end  and  clamped  end, 
and  the  16-crn  notches  have  larger  deviations  in  frequency  in  the  eccentrically  placed 
notches  rather  than  the  centered  notches.  The  difference  in  percent  deviations  is  a 
result  of  the  difference  in  notch  location  along  the  beam  thickness,  and  the  difference 
in  thickness. 

In  the  case  of  the  centrally  notched  beam,  the  modal  frequency  percent  de¬ 
viations  from  the  notchless  case  of  notches  with  equal  length  at  different  locations 
differ  less  than  5%  for  the  8-cm,  the  12-cm,  and  the  16-crn  notched  beams.  The 
largest  percent  deviation  occurs  in  the  4-cm  notched  beams  between  the  middle  and 
clamped-end  notched  beams,  which  differ  in  the  eighth  mode  by  7%.  Therefore,  the 
location  appears  more  difficult  to  diagnose  in  beams  containing  centered  notches. 
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CURVE-FIT  MODAL  FREQUENCIES 
Centrally  notched  beams 

Deviation  from  notchless  beam  (%) 

Mode 

4 

5 

6 

7 

8 

Notchless 

944 

1562 

2336 

3266 

4352 

Free  End 

4-cm 

909 

3.71% 

1497 

4.16% 

2233 

4.41% 

3117 

4.56% 

4149 

4.66% 

8-cm 

776 

17.80% 

1257 

19.53% 

1856 

20.55% 

2573 

21.22% 

3408 

21.69% 

12-cm 

711 

24.68% 

1150 

26.38% 

1695 

27.44% 

2346 

28.17% 

3103 

28.70% 

16-cm 

642.0 

31.99% 

1045 

33.10% 

1546 

33.82% 

2145 

34.32% 

2842 

34.70% 

Middle 

4-cm 

901 

4.56% 

1446 

7.43% 

2115 

9.46% 

2908 

10.96% 

3825 

12.11% 

8-cm 

742 

21.40% 

1223 

21.70% 

1828 

21.75% 

2557 

21.71% 

3410 

21.65% 

12-cm 

685 

27.44% 

1114 

28.68% 

1647 

29.49% 

2284 

30.07% 

3025 

30.49% 

16-cm 

584.0 

38.14% 

971 

37.84% 

1464 

37.33% 

2063 

36.83% 

2768 

36.40% 

Clamped  End 

4-cm 

882 

6.57% 

1441 

7.75% 

2138 

8.48% 

2973 

8.97% 

3946 

9.33% 

8-cm 

739 

21.72% 

1217 

22.09% 

1819 

22.13% 

2545 

22.08% 

3395 

21.99% 

12-cm 

677 

28.28% 

1098 

29.71% 

1621 

30.61% 

2246 

31.23% 

2973 

31.69% 

16-cm 

591.0 

37.39% 

976 

37.52% 

1461 

37.46% 

2046 

37.35% 

2731 

37.25% 

Table  4.6  Modal  frequencies  obtained  from  the  least-squares  curve-fit  of  the  ex¬ 
perimental  data  of  the  centrally  notched  beams. 
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Modal  Frequency  Percent  Deviation  of  4-cm  Centrally  Notched 

Beams 
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Figure  4.18  Modal  frequency  percent  deviation  of  4-cm  centrally  notched  beams. 


Modal  Frequency  Percent  Deviation  of  8-cm  Centrally  Notched 

Beams 
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Figure  4.19  Modal  frequency  percent  deviation  of  8-cm  centrally  notched  beams. 
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Modal  Frequency  Percent  Deviation  of  12-cm  Centrally  Notched 

Beams 
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Figure  4.20  Modal  frequency  percent  deviation  of  12-cm  centrally  notched  beams. 


Modal  Frequency  Percent  Deviation  of  16-cm  Centrally 
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Figure  4.21  Modal  frequency  percent  deviation  of  16-cm  centrally  notched  beams. 
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Mode 

5 

6 

7 

4-cm 

8-cm 

3%  to  8% 
9%  to  23% 

3%  to  10% 
10%  to  25% 

3%  to  10% 
10%  to  22% 

Table  4.7  Ranges  of  percent  deviation  from  notchless  case  for  4-cm  and  8-cm 
notched  beams  for  mode  5,  mode  6,  and  mode  7. 

Again,  observing  the  6th  mode,  if  the  deviation  is  between  4%  and  10%,  then 
the  notch  is  13%  of  the  beam.  If  the  deviation  is  between  20%  and  25%,  then  the 
notch  is  26%  of  the  beam.  If  the  deviation  is  between  26%  and  31%  of  the  beam, 
then  the  notch  is  39%  of  the  beam  length.  If  the  deviation  is  between  34%  and  38% 
then  the  notch  length  is  over  50%  of  the  beam. 

If  the  above  observation  is  considered  in  conjunction  with  the  observations  of 
the  eccentric  notches,  a  method  for  approximating  notch  size  despite  location  can  be 
determined.  A  frequency  deviation  of  3%  to  10%  indicates  a  notch,  which  is  around 
13%  of  the  beam  length.  Likewise,  a  deviation  of  10%  to  25%  indicates  a  notch 
length  around  26%  of  the  beam.  Longer  notches  cannot  be  distinguished  unless  the 
notch  location  is  known.  However,  the  evaluating  of  small  notches  is  more  critical  for 
practical  purposes.  The  other  modes  are  then  evaluated  and  the  results  are  displayed 
in  table  4.7 

4-4  Perel-Palazotto  Code 

The  Perel-Palazotto  code  was  first  run  for  the  notchless  beam  case.  The 
results  are  shown  compared  with  the  theoretical  natural  frequencies  shown  in  the 
following  table. 

The  theory  differs  by  less  than  3%.  The  Perel-Palazotto  code  was  run  for 
three  notched  beam  cases,  the  4-cm  and  8-cm  free-end  notched  beams  and  the  4-cm 
middle  notched  beam.  The  deviation  from  the  experimental  case  is  displayed  in  the 
following  table. 
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Mode 

Theory 

Perel-Palazotto 

Error 

1 

56.31 

57.44 

-2.01% 

2 

352.9 

359.1 

-1.77% 

3 

988.2 

1002 

-1.41% 

4 

1937 

1954 

-0.91% 

5 

3201 

3210 

-0.29% 

6 

4782 

4759 

0.47% 

7 

6678 

6578 

1.34% 

8 

8892 

8685 

2.32% 

Table  4.8  Comparison  of  theoretical  and  Perel-Palazotto  natural  frequencies  of  the 

notchless  beam. 


Modal  Frequencies 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

Exp. 

MATLAB 

A% 

52.92 

57.44 

7.87% 

330.9 

359.1 

7.87% 

1053 

1002 

-5.07% 

1869 

1954 

4.34% 

3093 

3210 

3.66% 

4590 

4759 

3.56% 

6306 

6589 

4.29% 

8450 

8685 

2.71% 

Free  End 

4-cm 

Exp. 

MATLAB 

A% 

52.81 

57.44 

8.06% 

328.1 

359.1 
8.63% 

1009 

1002 

-0.70% 

1828 

1954 

6.45% 

2972 

3210 

7.41% 

4322 

4758 

9.16% 

6178 

6586 

6.19% 

8134 

8678 

6.27% 

Free  End 

8-cm 

Exp. 

MATLAB 

A% 

52.50 

57.44 

8.60% 

325.0 

359.1 

9.50% 

987.5 

1002 

1.45% 

1731 

1953 

11.37% 

2903 

3196 

9.18% 

3966 

4684 

15.33% 

5481 

5857 

6.41% 

6881 

6717 

-2.44% 

Middle 

4-cm 

Exp. 

MATLAB 

A% 

51.88 

57.44 

9.68% 

331.3 

359.1 

7.74% 

812.5 

1002 

18.91% 

1853 

1954 

5.17% 

2887 

3209 

10.04% 

4309 

4757 

9.42% 

5844 

6585 

11.25% 

8159 

8671 

5.90% 

Table  4.9  Comparison  of  the  experimental  results  and  the  Perel-Palazotto  code 
results  for  selected  cases. 

The  results  of  the  Perel-Palazotto  code  agree  within  20%  of  the  experimental 
values.  In  general  the  experimental  natural  frequencies  suggest  a  lower  stiffness  than 
the  code  suggests.  This  difference  may  in  part  be  clue  to  a  small  difference  in  the 
material  properties  entered  in  the  code,  and  the  actual  material  properties. 
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4-5  Cobb-Durham  Code 


The  Cobb-Durham  code  was  first  run  for  the  notchless  case  and  compared 
with  the  theoretical  modal  frequencies,  which  can  be  seen  in  Table  4.10.  The  Cobb- 
Durham  code  differs  less  than  1%  from  the  theoretical  values. 


Mode 

Theory 

Cobb-Durham 

Error 

1 

56.31 

55.81 

0.88% 

2 

352.9 

349.7 

0.91% 

3 

988.2 

979.9 

0.90% 

4 

1937 

1919 

0.91% 

5 

3201 

3172 

0.90% 

6 

4782 

4739 

0.90% 

7 

6678 

6619 

0.89% 

8 

8892 

8813 

0.89% 

Table  4.10  Comparison  of  theoretical  and  Cobb-Durham  natural  frequencies  of  the 

notchless  beam. 


The  code  was  then  run  for  all  the  eccentric  notch  cases.  The  code  models 
the  section  containing  the  notch  as  a  beam  with  a  smaller  moment  of  inertia.  The 
computer  results  are  compared  with  the  experimental  results  in  Table  4.11.  The 
percent  difference  is  calculated  by 


%change 


(Experimental  modal  frequency)  —  (Computed  modal  frequency) 
(Experimental  modal  frequency) 


xl00% 

(4.4) 


The  Cobb-Durham  code  produces  close  results  for  the  4-cm  and  8-crn  notched 
beams.  The  4-cm  notched  beam  frequencies  are  predicted  within  5%  for  modes  higher 
than  mode  3.  The  modal  frequencies  of  the  8-crn  notched  beams  are  predicted  within 
10%  for  the  modes  higher  than  mode  3.  The  12-crn  notched  beam  frequencies  are 
also  within  10%  for  the  notch  located  in  the  middle.  The  clamped  end  and  free  end 
experimental  frequencies  are  much  lower  than  the  computed  frequencies,  as  well  as 
the  16-crn  notch  cases.  This  deviation  is  expected  since  the  code  does  not  account 
for  contact  between  crack  surfaces  and  allows  them  to  overlap  during  vibration. 
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EXPERIMENTAL  MODAL  FREQUENCIES 

Centrally  Located  Notch 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

Exp. 

MATLAB 

A% 

52.92 

55.81 

-5.46% 

330.9 

349.7 

-5.69% 

1053 

979.3 

7.00% 

1869 

1919 

-2.66% 

3093 

3172 

-2.57% 

4590 

4739 

-3.25% 

6306 

6619 

-4.96% 

8450 

8813 

-4.30% 

Free  End 

4-cm 

Exp. 

MATLAB 

A% 

52.81 

56.40 

-6.80% 

328.1 

351.1 
-7.01% 

1009 

967.0 

4.16% 

1828 

1840 

-0.66% 

2972 

2953 

0.64% 

4322 

4367 

-1.04% 

6178 

6140 

0.62% 

8134 

8256 

-1.50% 

8-cm 

Exp. 

MATLAB 

A% 

52.50 

56.69 

-7.98% 

325.0 

333.4 

-2.58% 

987.5 

840.7 

14.87% 

1731 

1664 

3.87% 

2903 

2809 

3.23% 

3966 

4098 

-3.33% 

5481 

5760 

-5.08% 

6881 

7762 

-12.80% 

12-cm 

Exp. 

MATLAB 

A% 

52.50 

56.25 

-7.14% 

306.3 

300.8 

1.80% 

942.5 

798.9 

15.24% 

1666 

1611 

3.30% 

2291 

2593 

-13.18% 

2550 

3975 

-55.88% 

3700 

5425 

-46.62% 

5023 

7367 

-46.67% 

16-cm 

Exp. 

MATLAB 

A% 

47.50 

54.13 

-13.96% 

295.0 

269.0 

8.81% 

352.5 

782.0 

-121.84% 

857.5 

1471 

-71.55% 

1485 

2472 

-66.46% 

2400 

3708 

-54.50% 

2805 

5100 

-81.82% 

3738 

6919 

-85.10% 

Middle 

4-cm 

Exp. 

MATLAB 

A% 

51.88 

53.48 

-3.08% 

331.3 

306.9 

7.36% 

812.5 

968.2 

-19.16% 

1853 

1749 

5.61% 

2887 

3067 

-6.25% 

4309 

4441 

-3.06% 

5844 

6264 

-7.19% 

8159 

8381 

-2.72% 

8-cm 

Exp. 

MATLAB 

A% 

50.94 

50.09 

1.67% 

328.1 

278.0 

15.27% 

775.0 

868.0 

-12.00% 

1878 

1720 

8.41% 

2459 

2691 

-9.43% 

4013 

4162 

-3.71% 

5744 

5848 

-1.81% 

7703 

7581 

1.58% 

12-cm 

Exp. 

MATLAB 

A% 

50.31 

47.78 

5.03% 

317.5 

269.6 
15.09% 

710.0 

786.0 

-10.70% 

1638 

1645 

-0.43% 

2678 

2658 

0.75% 

3850 

3839 

0.29% 

5473 

5471 

0.04% 

7308 

7413 

-1.44% 

16-cm 

Exp. 

MATLAB 

A% 

50.00 

45.40 

9.20% 

293.8 

264.1 

10.11% 

631.3 

733.1 

-16.13% 

906.3 

1502 

-65.73% 

985.0 

2567 

-160.61% 

1819 

3787 

-108.19% 

2865 

5151 

-79.79% 

3766 

6842 

-81.68% 

Clamped  End 

4-cm 

Exp. 

MATLAB 

A% 

50.94 

44.10 

13.43% 

318.2 

312.8 

1.70% 

955.7 

914.9 

4.27% 

1747 

1818 

-4.06% 

2858 

2998 

-4.90% 

4212 

4446 

-5.56% 

6279 

6193 

1.37% 

8169 

8274 

-1.29% 

8-cm 

Exp. 

MATLAB 

A% 

50.63 

38.59 

23.78% 

296.9 

309.4 

-4.21% 

871.9 

852.1 

2.27% 

1747 

1671 

4.35% 

2850 

2842 

0.28% 

3760 

4209 

-11.94% 

4031 

5790 

-43.64% 

5725 

7817 

-36.54% 

12-cm 

Exp. 

MATLAB 

A% 

49.38 

36.64 

25.80% 

267.5 

290.0 

-8.41% 

575.5 

805.7 

-40.00% 

880 

1644 

-86.82% 

1650 

2623 

-58.97% 

2593 

4013 

-54.76% 

3960 

5536 

-39.80% 

4975 

7385 

-48.44% 

16-cm 

Exp. 

MATLAB 

A% 

52.50 

35.57 

32.25% 

242.5 

256.2 

-5.65% 

352.5 

794.3 

-125.33% 

857.5 

1483 

-72.94% 

1485 

2530 

-70.37% 

2400 

3710 

-54.58% 

2805 

5216 

-85.95% 

3738 

6965 

-86.33% 

Table  4.11  Comparison  of  experimental  and  Cobb-Durham  results  for  the  eccen¬ 
trically  notched  beams. 


The  results  produced  for  the  centrally  notched  beams  are  displayed  in  Table 

4.12. 

The  natural  frequencies  of  the  centrally  notched  beams  show  even  better  agree¬ 
ment  between  the  code  and  the  experimental  results.  The  first  and  second  mode 
appear  to  show  some  large  deviation.  However,  the  code  modal  frequencies  of  mode 
4  and  higher  modes  are  within  10%  of  the  experimental  frequencies  for  most  cases. 

The  Cobb-Durham  code  does  not  account  for  the  nonlinearities  of  the  response, 

which  result  from  the  crack  surfaces  touching  when  the  sublaminates  are  vibrating  at 

different  frequencies.  The  nonlinear  nature  is  augmented  in  the  eccentrically  located 

notch,  which  is  the  reason  the  code  agreement  with  experimental  results  differs  by 

a  greater  amount  for  the  eccentrally  notched  beams  than  for  the  centrally  notched 
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beams.  However,  the  Cobb-Durham  code  provides  good  estimates  of  the  modal 
frequencies. 


Jh6  Notch  Depth 

The  comparison  of  the  centered  and  eccentric  notches  with  beams  of  the 
same  thickness  may  provide  more  information.  Using  the  Cobb-Durham  code,  modal 
frequencies  were  computed  for  beams  with  centered  notches  and  a  beam  thickness 
of  1/4- in.  These  frequencies  are  compared  with  the  eccentrically  notched  beams 
computed  frequencies  in  Figures  4.22  through  4.25. 

For  each  notch  length,  the  centrally  notched  beams  at  all  locations  have  higher 
modal  frequencies  than  all  three  of  the  eccentrically  notched  beams.  In  the  4-cm 


MODAL  FREQUENCIES 

Centrally  Located  Notch 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

Exp. 

MATLAB 

A% 

28.09 

27.90 

0.68% 

172.4 

174.9 

-1.43% 

484.1 

489.7 

-1.16% 

944.8 

959.5 

-1.56% 

1566 

1586 

-1.30% 

2342 

2370 

-1.20% 

3268 

3310 

-1.30% 

4359 

4407 

-1.09% 

Free  End 

4-cm 

Exp. 

MATLAB 

A% 

28.91 

28.54 

1.28% 

175.3 

175.9 

-0.36% 

471.7 

469.7 
0.43% 

931.8 

849.0 

8.93% 

1485 

1349 

9.16% 

2213 

2048 

7.47% 

3116 

2925 

6.13% 

4151 

3912 

5.76% 

8-cm 

Exp. 

MATLAB 

A% 

29.30 

28.74 

1.91% 

166.9 

152.5 

8.64% 

421.6 

358.1 

15.07% 

774.3 

757.0 

2.19% 

1360 

1220 

10.27% 

1779 

1798 

-1.05% 

2533 

2601 

-2.67% 

3477 

3353 

3.56% 

12-cm 

Exp. 

MATLAB 

A% 

29.04 

28.06 

3.37% 

154.3 

123.8 

19.77% 

406.3 

348.6 

14.19% 

658.2 

661.0 

-0.43% 

1090 

1124 

-3.09% 

1748 

1644 

5.97% 

2369 

2336 

1.39% 

3048 

3067 

-0.61% 

16- cm 

Exp. 

MATLAB 

A% 

28.26 

25.33 

10.37% 

146.8 

107.5 

26.76% 

336.6 

314.1 

6.68% 

561.5 

608.0 

-8.20% 

1046 

965 

7.78% 

1535 

1510 

1.63% 

2264 

2064 

8.82% 

2750 

2726 

0.87% 

Middle 

4-cm 

Exp. 

MATLAB 

A% 

27.73 

24.84 

10.42% 

170.6 

132.7 
22.20% 

456.2 

474.3 
-3.96% 

899.2 

819 

8.94% 

1536 

1449 

5.64% 

1978 

2137 

-8.06% 

3053 

2901 

4.98% 

3786 

4045 

-6.85% 

8-cm 

Exp. 

MATLAB 

A% 

28.13 

21.34 

24.14% 

169.0 

116.9 

30.83% 

390.9 

367.7 

5.93% 

725.8 

793 

-9.23% 

1283 

1191 

7.17% 

1847 

1776 

3.83% 

2477 

2637 

-6.46% 

3483 

3383 

2.88% 

12-cm 

Exp. 

MATLAB 

A% 

27.47 

19.42 

29.30% 

161.3 

111.9 

30.64% 

366.1 

314.5 

14.10% 

635.7 

689 

-8.43% 

1081 

1169 

-8.17% 

1729 

1636 

5.40% 

2378 

2240 

5.79% 

2991 

3118 

-4.25% 

16- cm 

Exp. 

MATLAB 

A% 

26.82 

17.72 

33.93% 

152.4 

106.3 

30.26% 

324.9 

289.5 

10.89% 

577.3 

584 

-1.23% 

1024 

1032 

-0.75% 

1403 

1582 

-12.73% 

1982 

2146 

-8.26% 

2820 

2757 

2.25% 

Clamped  End 

4-cm 

Exp. 

MATLAB 

A% 

27.73 

16.87 

39.17% 

169.9 

145.9 
14.13% 

474.9 

436.5 

8.08% 

965.5 

858 

11.12% 

1387 

1394 

-0.48% 

2077 

2075 

0.11% 

3093 

2934 

5.13% 

3946 

3970 

-0.60% 

8-cm 

Exp. 

MATLAB 

A% 

26.90 

13.71 

49.02% 

154.6 

137.6 
11.01% 

471.0 

367.0 

22.07% 

690.3 

762 

-10.35% 

1235 

1283 

-3.85% 

1763 

1830 

-3.76% 

2599 

2624 

-0.95% 

3404 

3524 

-3.53% 

12-cm 

Exp. 

MATLAB 

A% 

25.13 

12.75 

49.26% 

146.6 

116.4 

20.62% 

369.1 

352.7 

4.45% 

682.8 

694 

-1.68% 

1039 

1124 

-8.15% 

1656 

1722 

-3.96% 

2290 

2339 

-2.15% 

2949 

3206 

-8.71% 

16- cm 

Exp. 

MATLAB 

A% 

24.22 

12.25 

49.42% 

137.1 

94.30 

31.22% 

300.3 

325.4 
-8.36% 

563.7 

610.0 

-8.21% 

1026 

1001 

2.44% 

1449 

1549 

-6.90% 

1977 

2076 

-5.01% 

2756 

2841 

-3.10% 

Table  4.12  Comparison  of  experimental  and  Cobb-Durham  results  for  the  centrally 
notched  beams. 


4-28 


Computational  Modal  Frequency  Percent  Deviation  of  4-cnt 
Notched  1  4-inch  Beams 


Mode  llumhei 


Figure  4.22  Computational  modal  frequency  percent  deviation  of  4-cm  notched 
beams. 

case,  the  centrally  notched  beam  frequencies  are  at  least  2%  larger  than  eccentrically 
notched  beam  frequencies,  and  are  4%  larger  in  the  seventh  mode.  In  the  8-cm  case, 
the  centrally  notched  beam  frequencies  are  5%  to  7%  larger.  These  trends  indicate 


Computational  Modal  Fi  equency  Percent  Deviation  of  8-cm 


Eccentrically  Notched  1  4-inch  Beams 


Mode  llumhei 


Figure  4.23 


Computational  modal  frequency  percent  deviation  of  8-cm  notched 
beams. 
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the  possibility  of  estimating  the  notch  depth  with  experimentally  obtained  ranges  of 
the  percent  deviation  from  the  notchless  case. 


Figure  4.24 


Comp11t.1tion.1l  Mod.il  Fiequency  Percent  Deviation  of  12-cm 
Eccentrically  Notched  1. 4-inch  Beams 
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Computational  modal  frequency  percent  deviation  of  12-cm  notched 
beams. 
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Computational  Modal  Fiequency  Pei  cent  Deviation  of  16-cm 
Eccentrically  Notched  1  4-inch  Beams 
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Figure  4.25 


Computational  modal  frequency  percent  deviation  of  16-cm  notched 
beams. 
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4-7  Locations  of  maximum  shear 


The  notch  should  have  greater  impact  on  points  where  the  beam  is  sub¬ 
jected  to  greater  shear.  The  locations  of  maximum  shear  occur  at  nodes  of  modal 
curvature  inflection,  where  the  mode  shape  changes  from  concave  up  to  concave 
down,  or  vice  versa.  A  notch  that  intersects  one  of  these  locations  of  maximum 
shear  may  have  an  amplified  effect  on  the  modal  frequencies. 


In  Table  4.14,  the  change  in  frequency  as  the  notch  length  increases  in  length  by 
4  cm,  and  the  number  of  points  of  maximum  shear  crossed  by  the  notch  are  presented. 
The  large  percent  changes  should  correspond  to  a  difference  in  the  number  of  points 
crossed. 


%  Change 


(Frequency  of  longer  notched  beam)  —  (Frequency  of  shorter  notched  beam) 
(Frequency  of  shorter  notched  beam) 


(4.5) 


x  100% 


Several  modes  exhibit  the  expected  behavior.  In  the  beams  containing  notches 
at  the  clamped  end,  the  8-cm  notch  crosses  two  additional  nodes  in  the  seventh 
mode  and  one  additional  node  in  the  eighth  mode,  which  result  in  large  frequency 
changes  of  36%  and  30%  respectively.  The  12-crn  notch  crosses  an  additional  node 
in  mode  three  and  mode  5,  which  result  in  frequency  changes  of  34%  and  42%.  The 
16-cm  notch  crosses  an  additional  node  in  the  eighth  mode  and  two  additional  nodes 
in  the  seventh  mode,  which  result  in  25%  and  29%  changes.  However,  the  theory 
does  not  hold  up  for  some  instances.  The  8-cm  notch  crosses  additional  nodes  in  the 
fourth  and  fifth  modes,  but  the  modal  frequency  does  not  change.  The  8-cm  notch 
crosses  two  additional  modes  in  the  sixth  mode,  and  the  frequency  change  is  only 
11%.  In  the  case  of  the  12-crn  notch,  a  large  percentage  change  in  frequency  is  seen 
where  an  additional  node  is  not  crossed,  in  the  fourth  and  sixth  modes.  As  for  the 
16-cm  notch,  little  frequency  change  is  seen  where  additional  nodes  are  crossed  in 
the  fourth  and  sixth  modes. 


4-32 


EXPERIMENTAL  MODAL  FREQUENCIES 

Deviation  from  notchless  beam  (%) 

Number  of  nodes  crossed 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Notchless 

52.92 

330.9 

1053 

1869 

3093 

4590 

6306 

8450 

Free  End 

4-cm 

52.81 

328.1 

1009 

1828 

2972 

4322 

6178 

8134 

0.21% 

0.84% 

4.18% 

2.21% 

3.90% 

5.84% 

2.03% 

3.74% 

0 

0 

0 

0 

0 

0 

0 

0 

8-cm 

52.50 

325.0 

987.5 

1731 

2903 

3966 

5481 

6881 

0.79% 

1.77% 

6.22% 

7.40% 

6.14% 

13.59% 

13.08% 

18.57% 

0 

0 

0 

0 

1 

1 

1 

2 

12-cm 

52.50 

306.3 

942.5 

1666 

2291 

2550 

3700 

5023 

0.79% 

7.42% 

10.49% 

10.88% 

25.92% 

44.44% 

41.33% 

40.56% 

0 

0 

0 

1 

1 

2 

3 

3 

16-cm 

47.50 

295.0 

352.5 

857.5 

1485 

2400 

2805 

3738 

10.24% 

10.84% 

66.52% 

54.13% 

51.98% 

47.71% 

55.52% 

55.76% 

0 

0 

0 

1 

1 

1 

2 

2 

Middle 

4-cm 

51.88 

331.3 

812.5 

1853 

2887 

4309 

5844 

8159 

1.97% 

-0.13% 

1.96% 

0.87% 

6.66% 

6.12% 

7.33% 

3.44% 

0 

0 

0 

0 

1 

1 

1 

2 

8-cm 

50.94 

328.1 

775.0 

1878 

2459 

4013 

5744 

7703 

3.74% 

0.84% 

6.48% 

-0.46% 

20.49% 

12.57% 

8.92% 

8.84% 

0 

0 

1 

0 

1 

2 

2 

4 

12-cm 

50.31 

317.5 

710.0 

1638 

2678 

3850 

5473 

7308 

4.93% 

4.04% 

14.33% 

12.38% 

13.41% 

16.12% 

13.21% 

13.51% 

0 

0 

1 

1 

2 

2 

3 

4 

16-cm 

50.00 

293.8 

631.3 

906.3 

985 

1819 

2865 

3766 

5.52% 

11.20% 

23.82% 

51.52% 

68.15% 

60.37% 

54.57% 

55.43% 

0 

1 

2 

2 

3 

4 

5 

6 

Clamped  End 

4-cm 

50.94 

318.2 

955.7 

1747 

2858 

4212 

6279 

8169 

3.74% 

3.83% 

9.24% 

6.60% 

7.59% 

8.24% 

0.43% 

3.33% 

0 

0 

1 

1 

1 

1 

1 

2 

8-cm 

50.63 

296.9 

871.9 

1747 

2850 

3760 

4031 

5725 

4.33% 

10.27% 

17.20% 

6.54% 

7.85% 

18.08% 

36.08% 

32.25% 

0 

1 

1 

2 

2 

3 

3 

3 

12-cm 

49.38 

267.5 

575.5 

880 

1650 

2593 

3960 

4975 

6.69% 

19.15% 

45.35% 

52.92% 

46.65% 

43.51% 

37.21% 

41.12% 

0 

1 

2 

2 

3 

3 

3 

4 

16-cm 

52.50 

242.5 

352.5 

857.5 

1485 

2400 

2805 

3738 

0.79% 

26.71% 

66.52% 

54.13% 

51.98% 

47.71% 

55.52% 

55.76% 

0 

1 
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3 

3 

4 

5 

5 

Table  4.13  Modal  frequency  and  nodes  crossed  of  the  eccentrally  notched  beams. 


NODES  CROSSED  (Eccentric  Notch) 

(%)  Change 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Free  End 

4cm 

0 

0 

0 

0 

0 

0 

0 

0 

0.59% 

0.94% 

2.13% 

5.31% 

2.33% 

8.24% 

11.28% 

15.40% 

8cm 

0 

0 

0 

0 

1 

1 

1 

2 

0.00% 

5.75% 

4.56% 

3.76% 

21.07% 

35.70% 

32.50% 

27.00% 

12cm 

0 

0 

0 

1 

1 

2 

3 

3 

9.52% 

3.69% 

62.60% 

48.53% 

35.18% 

5.88% 

24.19% 

25.58% 

16cm 

0 

0 

0 

1 

1 

1 

2 

2 

Middle 

4cm 

0 

0 

0 

0 

1 

1 

1 

2 

1.81% 

0.97% 

4.62% 

-1.35% 

14.82% 

6.87% 

1.71% 

5.59% 

8cm 

0 

0 

1 

0 

1 

2 

2 

4 

1.24% 

3.23% 

8.39% 

12.78% 

-8.91% 

4.06% 

4.72% 

5.13% 

12cm 

0 

0 

1 

1 

2 

2 

3 

4 

0.62% 

7.46% 

11.08% 

44.67% 

63.22% 

52.75% 

47.65% 

48.47% 

16cm 

0 

1 

2 

2 

3 

4 

5 

6 

Clamped  End 

4cm 

0 

0 

1 

1 

1 

1 

1 

2 

0.61% 

6.69% 

8.77% 

0.00% 

0.28% 

10.73% 

35.80% 

29.92% 

8cm 

0 

1 

1 

2 

2 

3 

3 

3 

2.47% 

9.90% 

33.99% 

49.63% 

42.11% 

31.04% 

1.76% 

13.10% 

12cm 

0 

1 

2 

2 

3 

3 

3 

4 

-6.32% 

9.35% 

38.75% 

2.56% 

10.00% 

7.44% 

29.17% 

24.86% 

16cm 

0 

1 

2 

3 

3 

4 

5 

5 

Table  4.14  Percent  change  between  notch  size  and  number  points  of  maximum 
shear  crossed  in  the  eccentrally  notched  beam. 


In  the  case  of  the  mid-bar  notches,  similar  variation  is  found.  The  8-cm  notch 
crosses  additional  nodes  in  the  third,  sixth,  seventh,  and  eighth  modes,  where  the 
frequency  change  is  less  than  6%.  Furthermore,  the  largest  frequency  change  between 
the  4-cm  and  8-cm  notched  bar  is  seen  in  the  fifth  mode  where  no  additional  node 
is  crossed.  The  largest  frequency  change  from  the  8-cm  notch  to  the  12-cm  notch 
is  in  the  fourth  mode  where  an  additional  node  is  crossed.  An  additional  node  is 
also  crossed  in  the  seventh  mode,  but  the  frequency  change  is  less  than  5%.  The 
16-cm  notch  crosses  one  or  two  additional  notches  in  modes  four  through  eight  where 
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a  large  frequency  change  is  seen.  The  second  and  third  modes  only  show  a  small 
change,  even  though  an  additional  notch  is  crossed. 

The  beams  containing  notches  at  the  free  end  also  produce  similar  results.  In 
the  case  of  the  8-cm  notched  beam,  the  largest  changes  occur  when  additional  nodes 
are  crossed.  The  12-cm  notch  crosses  additional  nodes  in  the  sixth,  seventh,  and 
eighth  modes  where  a  large  change  is  evident.  However,  a  large  change  occurs  in 
the  fifth  mode  where  no  additional  node  is  crossed,  and  a  small  change  occurs  in  the 
fourth  mode  where  an  additional  node  is  crossed. 

Table  4.15  displays  the  number  of  nodes  crossed  by  the  centered  notch  and  the 
changes  in  modal  frequency  as  crack  size  increases  for  the  centrally  notched  beams. 

The  results  for  the  centrally  notched  beams  yield  the  same  ambiguous  corre¬ 
lation  and  are  described  in  detail  by  Reifsnyder.18 

When  the  centrally  and  eccentrically  notched  beams  are  compared,  the  eccen¬ 
trically  notched  beams  are  found  to  cross  either  as  many  or  more  nodes  than  the 
centrally  notched  beams  do  for  the  experimental  mode  shapes  for  the  4cm,  8cm, 
and  12cm  notched  beams.  The  eccentrically  notched  beams  cross  more  nodes  in  the 
fourth  and  eighth  modes  for  the  middle  4cm  notch  and  in  the  third,  seventh,  and 
eighth  modes  for  the  clamped  end  4cm  notch.  The  free  end  4cm  notch  crosses  the 
same  number  of  nodes  in  both  cases.  The  eccentrically  notched  beams  cross  more 
nodes  in  the  fifth,  sixth,  and  eighth  modes  for  the  free  end  8cm  notch,  in  the  fourth, 
sixth,  and  eighth  modes  for  the  middle  8cm  notch,  and  in  the  fourth  through  eighth 
modes  for  the  clamped  end  8cm  notch.  Finally,  the  eccentrically  notched  beams 
cross  more  nodes  in  the  fourth,  sixth,  seventh,  and  eighth  modes  for  the  free  end 
12cm  notch,  in  the  eighth  mode  for  the  middle  12cm  notch,  and  in  the  third,  fourth, 
fifth,  seventh,  and  eighth  modes  for  the  clamped  end  12cm  notch. 

Thus  the  eighth  mode  appears  to  lend  some  indication  of  notch  depth  through 
the  modal  shape.  If  the  size  and  location  is  known,  then  the  depth  (if  limited  to  two 
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NODES  CROSSED  (Centered  Notch) 

(%)  Change 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Free  End 

4cm 

0 

0 

0 

0 

0 

0 

0 

0 

-1.35% 

4.76% 

10.62% 

16.90% 

8.44% 

19.61% 

18.70% 

16.25% 

8cm 

0 

0 

0 

0 

0 

0 

1 

1 

0.89% 

7.57% 

3.64% 

15.00% 

19.81% 

1.74% 

6.49% 

12.32% 

12cm 

0 

0 

0 

0 

1 

1 

2 

2 

2.69% 

4.88% 

17.15% 

14.70% 

4.04% 

12.20% 

4.45% 

9.79% 

16cm 

0 

0 

0 

1 

1 

2 

2 

3 

Middle 

4cm 

0 

0 

0 

1 

1 

1 

1 

1 

-1.44% 

0.92% 

14.33% 

19.29% 

16.45% 

6.62% 

18.87% 

7.99% 

8cm 

0 

0 

1 

1 

1 

1 

2 

2 

2.35% 

4.54% 

6.33% 

12.41% 

15.77% 

6.35% 

4.01% 

14.13% 

12cm 

0 

0 

1 

1 

2 

2 

3 

3 

2.37% 

5.52% 

11.27% 

9.19% 

5.21% 

18.85% 

16.63% 

5.71% 

16cm 

0 

0 

1 

2 

2 

3 

3 

4 

Clamped  End 

4cm 

0 

0 

0 

1 

1 

1 

0 

1 

3.00% 

8.99% 

0.82% 

28.50% 

10.95% 

15.12% 

15.95% 

13.75% 

8cm 

0 

1 

1 

1 

1 

2 

1 

2 

6.57% 

5.17% 

21.62% 

1.09% 

15.88% 

6.06% 

11.91% 

13.36% 

12cm 

0 

1 

1 

1 

2 

3 

2 

3 

3.62% 

6.50% 

18.65% 

17.44% 

1.28% 

12.52% 

13.66% 

6.56% 

16cm 

0 

1 

1 

2 

3 

3 

3 

4 

Table  4.15  Percent  change  between  notch  size  and  number  points  of  maximum 
shear  crossed  in  the  centrally  notched  beam. 


possibilities)  can  be  determined  by  the  following.  Considering  the  free  end,  if  an  8 
cm  notch  crosses  1  node,  it  is  centrally  located,  and  if  it  crosses  2  nodes,  the  notch 
is  eccentrically  located.  Two  nodes  crossed  by  a  12cm  notch  indicates  a  centered 
notch,  and  three  nodes  crossed  indicate  an  eccentric  notch.  Considering  the  middle 
and  clamped  end  notches,  if  a  4  cm  notch  crosses  1  node,  it  is  centrally  located, 
whereas  if  it  crosses  2  nodes  it  is  eccentrically  located.  Two  nodes  crossed  by  an  8 
cm  notch  indicate  a  centered  notch,  whereas  three  nodes  crossed  indicate  an  eccentric 
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notch.  Likewise,  if  a  12-cm  notch  crosses  three  nodes  it  is  centered,  and  if  a  12-cm 
notch  crosses  four  nodes  it  is  eccentric. 


Jh8  Locations  of  Maximum  Curvature 


Regions  of  maximum  curvature,  or  anti-nodes,  are  also  analyzed.  Curva¬ 


ture  at  a  point  is 


// 


V 


M 

EC 


(4.6) 


where  M  is  the  bending  moment,  E  is  the  modulus  of  elasticity,  and  /  is  the  second 
moment  of  inertial  of  the  cross-sectional  area.  Damage  results  in  a  decrease  in  the 
second  moment  of  inertia,  which  should  then  result  in  an  increase  in  curvature.  The 
table  below  shows  the  number  of  anti-nodes  crossed  by  the  notch  at  each  mode  for 
each  case. 


The  clamped-end  notched  beams  demonstrate  the  following  behavior.  The 
greatest  changes  in  frequency  from  the  4-cm  notched  beam  to  the  8-cm  notched 
beam  correspond  to  the  modes  where  additional  anti-nodes  are  crossed,  namely  the 
third,  sixth,  seventh,  and  eighth  modes.  The  12-cm  notch  crosses  additional  anti¬ 
nodes  in  the  third,  fifth,  and  sixth  modes,  where  large  changes  in  frequency  are 
demonstrated.  However,  the  frequency  change  is  large  in  the  fourth  mode  where  no 
additional  anti-node  is  crossed,  and  is  less  than  15%  for  the  eighth  and  second  modes 
where  an  additional  anti-node  is  crossed.  The  16-crn  notch  crosses  one  additional 
anti-node  in  the  eight  mode  and  two  additional  anti-nodes  in  the  seventh  mode 
where  the  frequency  changes  are  25%  and  30%.  Nonetheless,  a  large  frequency 
change  is  observed  where  no  additional  anti-node  is  crossed  in  the  third  mode,  and 
the  frequency  change  is  insignificant  where  an  additional  anti-node  is  crossed  in  the 
fifth  mode. 


In  the  mid-located  notched  beam,  the  largest  frequency  changes  from  the  4-cm 
notched  beam  to  the  8-cm  notched  beam  occur  when  the  notch  crosses  an  additional 
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ANTI-NODES  CROSSED  (Eccentric  Notch) 

(%)  Change 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Free  End 

4cm 

0 

0 

0 

0 

0 

1 

0 

1 

0.59% 

0.94% 

2.13% 

5.31% 

2.33% 

8.24% 

11.28% 

15.40% 

8cm 

0 

0 

0 

1 

1 

2 

2 

2 

0.00% 

5.75% 

4.56% 

3.76% 

21.07% 

35.70% 

32.50% 

27.00% 

12cm 

0 

0 

1 

2 

2 

3 

3 

4 

9.52% 

3.69% 

62.60% 

48.53% 

35.18% 

5.88% 

24.19% 

25.58% 

16cm 

0 

1 

1 

1 

1 

2 

2 

3 

Middle 

4cm 

0 

1 

0 

1 

1 

1 

1 

1 

1.81% 

0.97% 

4.62% 

-1.35% 

14.82% 

6.87% 

1.71% 

5.59% 

8cm 

0 

1 

1 

1 

2 

2 

3 

2 

1.24% 

3.23% 

8.39% 

12.78% 

-8.91% 

4.06% 

4.72% 

5.13% 

12cm 

0 

1 

2 

1 

2 

3 

2 

3 

0.62% 

7.46% 

11.08% 

44.67% 

63.22% 

52.75% 

47.65% 

48.47% 

16cm 

0 

1 

2 

2 

3 

3 

4 

3 

Clamped  End 

4cm 

0 

0 

0 

0 

1 

1 

2 

1 

0.61% 

6.69% 

8.77% 

0.00% 

0.28% 

10.73% 

35.80% 

29.92% 

8cm 

0 

0 

1 

2 

1 

2 

3 

3 

2.47% 

9.90% 

33.99% 

49.63% 

42.11% 

31.04% 

1.76% 

13.10% 

12cm 

0 

1 

2 

2 

2 

3 

2 

4 

-6.32% 

9.35% 

38.75% 

2.56% 

10.00% 

7.44% 

29.17% 

24.86% 

16cm 

0 

1 

2 

2 

3 

3 

4 

5 

Table  4.16  Percent  change  between  notch  size  and  number  anti- nodes  crossed  in 
the  eccentrically  notched  beam. 


anti-node,  although  another  anti-node  is  crossed  in  the  seventh  mode  where  negligible 
frequency  change  is  seen.  The  12-cm  notch  crosses  additional  anti-nodes  in  modes 
five  through  eight,  and  the  corresponding  frequency  changes  are  large.  The  behavior 
of  the  third  and  fourth  modes  do  not  follow  this  trend. 

The  following  table  displays  the  number  of  nodes  crossed  by  the  centered  notch 
and  the  changes  in  modal  frequency  as  crack  size  increases  for  the  centrally  notched 
beams18 
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ANTI-NODES  CROSSED  (Centered  Notch) 

(%)  Change 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

Free  End 

4cm 

0 

0 

0 

0 

0 

0 

0 

0 

-1.35% 

4.76% 

10.62% 

16.90% 

8.44% 

19.61% 

18.70% 

16.25% 

8cm 

0 

0 

0 

0 

1 

1 

1 

1 

0.89% 

7.57% 

3.64% 

15.00% 

19.81% 

1.74% 

6.49% 

12.32% 

12cm 

0 

0 

0 

1 

1 

2 

2 

2 

2.69% 

4.88% 

17.15% 

14.70% 

4.04% 

12.20% 

4.45% 

9.79% 

16cm 

0 

0 

1 

1 

2 

2 

3 

3 

Middle 

4cm 

0 

0 

0 

1 

1 

1 

1 

1 

-1.44% 

0.92% 

14.33% 

19.29% 

16.45% 

6.62% 

18.87% 

7.99% 

8cm 

0 

0 

1 

1 

1 

2 

2 

2 

2.35% 

4.54% 

6.33% 

12.41% 

15.77% 

6.35% 

4.01% 

14.13% 

12cm 

0 

1 

0 

1 

2 

2 

3 

3 

2.37% 

5.52% 

11.27% 

9.19% 

5.21% 

18.85% 

16.63% 

5.71% 

16cm 

0 

1 

1 

2 

2 

3 

3 

4 

Clamped  End 

4cm 

0 

0 

1 

0 

0 

1 

1 

2 

3.00% 

8.99% 

0.82% 

28.50% 

10.95% 

15.12% 

15.95% 

13.75% 

8cm 

0 

1 

1 

1 

1 

2 

2 

3 

6.57% 

5.17% 

21.62% 

1.09% 

15.88% 

6.06% 

11.91% 

13.36% 

12cm 

0 

1 

1 

1 

2 

2 

3 

3 

3.62% 

6.50% 

18.65% 

17.44% 

1.28% 

12.52% 

13.66% 

6.56% 

16cm 

0 

1 

2 

1 

2 

3 

4 

4 

Table  4.17  Percent  change  between  notch  size  and  number  anti-nodes  crossed  in 
the  centrally  notched  beam. 


The  location  of  maximum  curvature  in  the  eighth  mode  reveals  information 
about  the  depth  of  a  notch  located  at  the  free  end.  If  a  4  cm  notch  crosses  1  anti- 
node,  it  is  eccentric,  and  if  it  crosses  no  anti-nodes,  it  is  centered.  Two  anti-nodes 
crossed  by  an  8  cm  notch  indicate  an  eccentric  notch,  whereas  one  anti-node  indicates 
a  centered  notch.  Likewise,  if  a  12  cm  notch  crosses  four  anti-nodes,  the  notch  is 
eccentric,  whereas  if  it  crosses  two,  it  is  centered. 


4-39 


4.9  ABAQUS  Model 

Two  centrally  notched  beam  specimens  were  modeled  using  the  ABAQUS 
CAE  package,  the  beams  containing  a  4  cm  notch  and  an  8  cm  notch  at  the  clamped 
end.  The  ABAQUS  approach  requires  that  each  case  be  individually  modeled  and 
meshed,  in  contrast  to  the  MATLAB  programs.  The  4  cm  notch  was  modeled  by 
Capt  Aaron  Reifsnyder.  The  8  cm  notch  was  modeled  by  Mr  Vikram  Dhruva.  Ap¬ 
proximately  14,000  8-noded  linear  hexahedral  brick  elements  were  used  to  model  the 
4  cm  notched  beam.  Approximately  60%  of  the  total  elements  were  concentrated 
in  the  notched  region.  A  mesh  of  66,600  elements  was  constructed  for  the  8  cm 
notched  beam,  also  concentrating  elements  in  the  notched  region.  The  beam  charac¬ 
teristics  required  for  ABAQUS  were  a  modulus  of  elasticity  of  10  x  106  psi,  a  density 
of  2.666  x  10-4  slugs  per  cubic  feet,  and  a  poisson’s  ratio  of  0.3.  The  piezoelectric 
patch  was  ignored  for  both  cases.  ABAQUS  computes  bending,  torsional,  and  mixed 
modes,  so  the  mode  shapes  were  once  again  used  to  distinguish  the  pure  bending 
modes.  The  bending  modes  are  displayed  in  the  table  below. 


MODAL  FREQUENCIES 

Mode 

1 

2 

3 

4 

5 

6 

7 

8 

4  cm 

Clamped  End 

Exp. 

ABAQUS 

A% 

28.65 

28.6 

-017% 

174.7 

175.5 

0.46% 

489.5 

498.2 

1.78% 

995.4 

995 

-0.04% 

1439 

1502 

4.38% 

2137 

2151 

0.66% 

3191 

3189 

-0.06% 

4105 

4311 

5.02% 

8  cm 

Clamped  End 

Exp. 

ABAQUS 

A% 

27.73 

26.86 

-3.13% 

159.4 

159.5 
0.08% 

485.5 

463.7 

-4.49% 

711.7 

713.8 
0.30% 

1274 

1309 

2.80% 

1818 

1842 

1.33% 

2680 

2646 

-1.26% 

3509 

3445 

-1.84% 

Table  4.18  Comparison  of  experimental  ABAQUS  modal  frequencies  for  the  4-cm 
and  8-cm  clamped-end  notched  beams. 


The  ABAQUS  results  are  within  5%  of  the  experimental  results.  This  deviation 
could  be  further  reduced  by  a  mesh  convergence  study.  ABAQUS  could  be  useful  in 
studying  the  modal  characteristics  of  the  eccentrally-located  notched  beam,  where 
the  nonlinear  behavior  is  augmented.  The  first  eight  modal  shapes  computed  by 
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the  ABAQUS  method  for  the  centrally  located  8-cm  clamped-end  notched  beam  are 
shown  below. 


.vi» 
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Figure  4.26  ABAQUS  bending  mode  1. 
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Figure  4.27  ABAQUS  bending  mode  2. 
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Figure  4.28  ABAQUS  bending  mode  3. 
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Figure  4.29  ABAQUS  bending  mode  4. 
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Figure  4.30  ABAQUS  bending  mode  5. 


Figure  4.31  ABAQUS  bending  mode  6. 
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Figure  4.32  ABAQUS  bending  mode  7. 


Figure  4.33  ABAQUS  bending  mode  8. 
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V.  Conclusions  and  Recommendations  for  Further  Study 

The  present  study  investigated  the  natural  frequency  response  of  A1  2024-T3 
cantielever  beams  containing  notches  of  varying  length  and  location.  The  diagnostic 
approaches  used  were  modal  analysis  and  finite  element  methods.  The  natural  fre¬ 
quency  trends  were  examined  according  to  notch  length,  location  along  the  length 
of  the  beam,  and  location  with  respect  to  the  beam  axis.  The  following  conclusions 
have  been  drawn. 

1.  The  natural  frequency  of  the  beam  decreased  as  the  notch  length  increases. 

2.  The  percent  decrease  from  the  notchless  case  of  the  natural  frequency  of 
beams  with  inclusions  was  generally  greater  for  a  beam  containing  notch  located  on 
the  beam  axis  than  off-axis. 

3.  Eccentric  notches  and  centered  notches  can  be  sized  using  the  natural 
frequency  deviation  from  the  notchless  case. 

4.  The  trend  of  the  natural  frequency  deviation  from  the  notchless  case  pro¬ 
vides  information  on  the  location  of  the  notch  in  the  case  of  eccentrically  located 
notches. 

6.  The  eccentric  notches  intersect  more  points  of  maximum  shear  and  points 
of  maximum  curvature  of  the  modal  shapes  than  the  centered  nodes. 

7.  The  acoustic  horn  introduced  interference  in  the  FRF  plot  of  the  notch¬ 
less  case  that  was  not  present  in  the  FRF  plot  of  the  piezoelectric-actuator  excited 
notchless  beam.  However,  it  is  difficult  to  attribute  the  interference  specifically  to 
the  acoustic  horn  or  to  notch  contact. 

8.  The  Cobb-Durham  code  is  a  simple  finite  element  model,  which  produces 
good  estimates  for  the  natural  frequencies  of  the  notched  beams  for  both  the  eccentric 
and  the  centered  notch  cases. 
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9.  ABAQUS  is  an  effective  modeling  tool  to  model  the  frequency  response  of 
beams  containing  notches. 

Further  experimentation  should  be  conducted  to  characterize  the  behavior  of 
the  acoustic  horn.  The  acoustic  horn  is  a  less  expensive  tool  than  the  piezoelectric 
actuator,  but  it  can  cause  interference  in  sensitive  measurements.  If  this  influence 
were  characterized,  the  horn  would  be  available  for  more  applications  and  would 
provide  better  results. 

Likewise,  the  notch  sublaminates  need  to  be  modeled,  to  determine  that  notch 
surface  contact  does  indeed  occur,  and  to  determine  the  circumstances  under  which 
interference  may  be  attributed  to  this  surfact  contact. 

With  respect  to  this  study,  the  experiment  may  be  performed  using  centrally 
notched  beams  with  the  same  thickness  as  the  eccentric  notches  in  this  study.  The 
direct  comparison  would  allow  verification  of  the  code  predictions,  as  well  as  perhaps 
more  insight  into  the  influence  of  notch  depth. 

Other  methods  of  locating  notches  are  possible.  An  experiment  could  be  con¬ 
ducted  to  measure  the  frequency  response  in  the  time  domain  and  use  a  spectral 
element  model  to  observe  the  stress  wave  as  it  propagates  down  the  length  of  the 
beam.  The  stress  wave  is  magnified  at  the  crack  tips,  and  thereby  can  indicate  both 
notch  size  and  location.  To  reexamine  the  notch  depth,  the  beams  of  this  experi¬ 
ment  could  be  tested  with  the  lower  sublaminate  facing  the  laser  vibrometer.  This 
orientation  may  results  in  different  measured  modal  frequencies  and  thereby  may 
offer  a  model  for  notch  depth. 
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Appendix  A.  Mode  Shapes  and  Frequencies  of  the  Eccentrically 

Notched  Beams 


Figure  A.l 
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First  eight  mode  shapes  and  modal  frequencies  of  the  notchless  beam. 
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Figure  A. 2 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  4-cm  notch  at  the  clamped  end. 
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Figure  A. 3 
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First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  8-cm  notch  at  the  clamped  end. 
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Figure  A. 4 
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First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  12-crn  notch  at  the  clamped  end. 


A-4 


„  iq-4  Mode152.50  Hz  x  10?  Mo<te2  242.5  Hz 


i 

0 

-1 

1 

x  10*  Mode3  352.5  Hz 

*  to'5  Mode4  057  5  Hz 

V 

0 

•1 

2 

□ 

•2 

5 

0 

-5 

a y  ■ 

„  IQ*  Mode5  1405  Hz 

*  io',u  Mode6  2400  Hz 

AA 

WA 

„  iQ*"  Mode7  2005  Hz 

x  10  ’u  ModeB  3738  Hz 

m 

w 

Figure  A. 5 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  16-crn  notch  at  the  clamped  end. 


A-5 


x  i0'J  Model  5188  Hz 


x  10"*  M«d«2  331.3  Hz 


— - — 

2 

a 

-2 

5 

0 

.c 

Mode3  812.5  Hz 

x  io-*  Mode4  18G3  Hz 

^\7 

x  *  Mode5  2688  Hz 

5 

x  iq  lu  Mode6  4309  Hz 

A/V 

0 

-5 

C 

yw 

x10"  Mode/ 5844  Hz 

x  ,0U  Mode8  0159  Hz 

0 

yxyxyvy 

-5 

Figure  A. 6 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  4-cm  notch  in  the  middle. 
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Figure  A.  7 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  8-cm  notch  in  the  middle. 
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Figure  A. 8 


„  iq-4  Model  50.31  Hz 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  12-cm  notch  in  the  middle. 
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Figure  A. 9 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  16-cm  notch  in  the  middle. 
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Figure  A. 10 
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First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  4-cm  notch  at  the  free  end. 


A-10 


X  in'®  Model  52  50  Hz 


'1x10^  Mode5  2903  Hz 


x  io'7  Mod«2  325.0  Hz 


Figure  A. 11 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  8-cm  notch  at  the  free  end. 
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Figure  A. 12 
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First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  12-cm  notch  at  the  free  end. 
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Figure  A. 13 


First  eight  mode  shapes  and  modal  frequencies  of  the  beam  containing 
the  16-cm  notch  at  the  free  end. 
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Appendix  B.  Cobb-Durham  Code 

7.  FEM  analytical  solution 
°/„  Written  by  Dr  Richard  Cobb 
7„  Modified  by  Frances  Durham 

7„  This  code  computes  the  natural  frequencies  and  mode  shapes  of  a 
°/„  cantilevered  beam  containing  an  eccentric  notch. 

beam_length  =  12/12;  °/„  ft 
h  =  .25/12  ;  %  ft 

h_u  =  h/4-0. 006/12;  h_l  =  3*h/4-0 . 006/12 ; 
b  =  1.5/12;  %  ft 

E  =  10000000*12*12 ;  %  lb-f/ft“2  (10000  ksi) 
rho  =  0 . 0975/32 . 174*12~3 ;  %  slug/ft “3 


n_els  =  36;  n_nodes  =  n_els  +  1;  j  =  n_nodes*2; 
c  =  .000001;  */.  assume  lightly  damped 

A(l:n_els)  =  b*h;  I(l:n_els)  =  b*(h~3)/12;  L(l:n_els)  = 
beam_length/n_els ; 

n_crack_tip_els  =  3;  crack_length_cm  =  input ( 5 Length  of  crack  in  cm 

’ );  crack_length  =  crack_length_cm/2. 54/12;  location  = 

input ( ’Location  of  crack  (1)  clamped  end  (2)  middle  (3)  free  end: 

’ );  n_crack_els  =  crack_length_cm*6/4;  L_crack  = 

crack_length/n_crack_els ;  L_crack_tip_els=L_crack;  if  (location==l) 
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L_undamaged  =  (beam_length  -  crack_length-n_crack_tip_els*L_crack_tip_els)/(n_el£ 
L ( 1 :n_crack_els)  =  L_crack; 

L(n_crack_els+1 :n_crack_els+n_crack_tip_els)  =  L_crack_tip_els ; 
L(n_crack_els+n_crack_tip_els+l :n_els)  =  L_undamaged; 
crackell=l ; 
crackel2=n_crack_els ; 
for  k=l :n_crack_els 

I (k)  =  b* (h_u~3+h_l~3) / 1 2 ; 

A(k)  =  b*(h_u+h_l); 

end 

slope  =  (I (n_crack_els+l)-I (n_crack_els) )/(n_crack_tip_els+l) ; 
for  k=l :n_crack_tip_els 

I (n_crack_els+k)  =  I(n_crack_els)+k*slope; 

end 

end  if  (location==2) 

L_undamaged  =  (beam_length-crack_length-2*n_crack_tip_els*L_crack_tip_els) / (n_el£ 
n_els_undamaged  =  n_els-n_crack_els-2*n_crack_tip_els ; 

L ( 1 : n_els_undamaged/2) =L_undamaged ; 

L(n_els_undamaged/2+l :n_els_undamaged/2+n_crack_tip_els)  =  L_crack_tip_els ; 
L(n_els_undamaged/2+n_crack_tip_els+l : n_els_undamaged/2+n_crack_tip_els+n_crack_€ 
L(n_els-n_els_undamaged/2-n_crack_tip_els+l :n_els-n_els_undamaged/2)  =  L_crack_tj 
L(n_els-n_els_undamaged/2 :n_els)  =  L_undamaged ; 
crackell=n_els_undamaged/2+n_crack_tip_els+l ; 
crackel2=crackell+n_crack_els ; 

for  k=n_els_undamaged/2+n_crack_tip_els+l :n_els_undamaged/2+n_crack_tip_els+n_cr£ 
I (k)  =  b* (h_u~3+h_l~3) / 12 ; 

A(k)  =  b* (h_u+h_l) ; 

end 


B-2 


slope  =  (I(l)-I(n_els_undamaged/2+n_crack_tip_els+l))/(n_crack_tip_els+2) ; 
for  k=l :n_crack_tip_els 

I (n_els_undamaged/2+k)  =  I(n_els_undamaged/2)-k*slope; 

I (n_els_undaraaged/2+n_crack_tip_els+k)  =  I (n_els_undaraaged/2+n_crack_tip_els^ 

end 

end  if  (location==3) 

L_undamaged  =  (beam_length  -  crack_length- (n_crack_tip_els+l) *L_crack_tip_els) / (r 
L ( 1 :n_els-n_crack_tip_els-n_crack_els-2)=L_undamaged; 

L(n_els-n_crack_tip_els-n_crack_els-l :n_els-n_crack_els-2)=L_crack_tip_els ; 
L(n_els-n_crack_els-l :n_els-l)=L_crack; 

L(n_els)=L_crack_tip_els ; 
crackel2=n_els-l ; 
crackell=crackel2-n_crack_els ; 
for  k=n_els-n_crack_els-l :n_els-l 
I (k)  =  b* (h_u~3+h_l~3) / 12 ; 

A(k)  =  b* (h_u+h_l) ; 

end 

slope  =  (I(l)-I(n_els-l))/(n_crack_tip_els+l) ; 
for  k=l :n_crack_tip_els 

I (n_els-n_crack_els-l-k)=I (n_els-l)+(n_crack_tip_els-k) *slope ; 

end 

end  pos(l)=0;  for  k=l:n_els-l 
pos(k+l)=pos(k)+L(k) ; 

end 

M  =  zeros(j);K  =  zeros(j); 

for  i=l:2:j-3 
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k=(i+l)/2; 

mel  =  (rho*A(k)*L(k)/420)* [156  22*L(k)  54  -13*L(k) ; . . . 

22*L(k)  4*(L(k)~2)  13*L(k)  -3* (L(k) ~2) ; . . . 
54  13*L(k)  156  -22*L(k) ; . . . 

-13*L(k)  -3*(L(k)“2)  -22*L(k)  4*(L(k)~2)] 

kel  =  (E*I (k)/ (L(k) ~3) ) *  [12  6*L(k)  -12  6*L(k);... 

6*L(k)  4*(L(k)“2)  -6*L(k)  2* (L(k) ~2) ; . . . 

-12  -6*L(k)  12  -6*L(k) ; . . . 

6*L(k)  2* (L(k) ~2)  -6*L(k)  4*(L(k)“2)]; 

M(i : i+3, i : i+3)  =  M(i : i+3 , i : i+3)  +  mel; 

K(i : i+3, i : i+3)  =  K(i : i+3 , i : i+3)  +  kel; 

end 


M  =  M(3: j ,3: j) ;  %  Apply  BC 

K=  K(3:j,3:j);  C  =  K*c;  F=zeros (n_els*2 , 1) ; 

F (6) =1 ;  %  patched  creates  moment  at  1.5  in  ->  node  3  ->  DOF  6 


°/  V  7  °/  7  °/  °/  7  7  °/  7 

/o  /o  /o  /  o  /o  /o  /o  /o  /o  /o  /o 

A=  [zeros (n_els*2 ,n_els*2)  eye (n_els*2) ; -M~-1*K  -M“-1*C]  ; 
[asize , other] =size (A) ;  C=eye(n_els*4) ;  C=C(1 :n_els*2, : ) ; 
[V,D] =eig(A) ;  for  i=l: asize 
Doh(i)=D(i,i) ; 

Dohabs (i)=abs (Doh(i) ) /2/pi ; 


end 

°/o 
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figured)  subplot  (4, 2  >  1)  M1=V( :  ,  asize) 

mlhz=abs (D (asize , asize) ) /2/pi ;  model=Ml (1 : 2 : 2*n_els) ; 

model=real (model/raax(model) ) ; 

plot (pos , model) ; title (sprintf ( ’Model  %3.3f  ’ ,mlhz)) 
set(gca,  ’XTick’  ,  [] ) 

°/o 

subplot (4, 2 ,2)  M2=V( : ,asize-2)  ra2hz=abs(D(asize-2,asize-2))/2/pi 
mode2=M2 (1 : 2 : 2*n_els) ;  raode2=real (mode2/max(mode2) ) 
plot (pos ,mode2) ; title (sprintf ( ’Mode2  703.3f  ’  ,m2hz')') 
set(gca,  ’XTick’  ,  [] ) 

i=0;  p=0;  mode2_crossedanti=0 ;  raode2_crossednode=0 ;  for  k=2:n_els-l 
if  mode2(k-l)<=mode2(k)  &  mode2(k+l)<=mode2(k) 
i=i+l ; 

mode2_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

mode2_crossedanti=mode2_crossedanti+l ; 

end 

elseif  mode2(k-l)>=mode2(k)  &  mode2(k+l)>=mode2(k) 
i=i+l ; 

mode2_maxcurve (i)=pos (k) ; 
if  k>=crackell  &  k<=crackel2 

mode2_crossedanti=mode2_crossedanti+l ; 

end 

end 

if  mode2(k)<=0  &  mode2(k+l)>=0 
p=p+l; 
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mode2_nodes (p) =pos (k) ; 
mode2_crossednode=mode2_crossednode+l ; 
elseif  mode2(k)>=0  &  mode2(k+l)<=0 
p=p+l; 

mode2_nodes(p)=pos(k) ; 
mode2_crossednode=mode2_crossednode+l ; 

end 

end 


°/o 

subplot (4, 2, 3)  M3=V( : , asize-4) ;  m3hz=abs (D (asize-4 , asize-4) ) /2/pi 
mode3=M3(l : 2 : 2*n_els) ;  raode3=real (mode3/max(mode3) ) ; 
plot (pos ,mode3) ; title (sprintf ( ’Mode3  703.2f  ’  ,m3hz)) 
set(gca,  ’XTick’  ,  [] ) 

°/0set  (gca,  ’YTick'  ,  [] ) 

1=0 ;  p=0;  mode3_crossedanti=0 ;  mode3_crossednode=0 ;  for  k=3:n_els-l 
if  mode3(k-l)<=mode3(k)  &  mode3(k+l)<=mode3(k) 
i=i+l ; 

mode3_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

mode3_crossedanti=mode3_crossedanti+l ; 

end 

elseif  mode3(k-l)>=mode3(k)  &  mode3(k+l)>=mode3(k) 
i=i+l ; 

if  k>=crackell  &  k<=crackel2 

raode3_crossedanti=mode3_crossedanti+l ; 

end 
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mode3_maxcurve (i)=pos (k) ; 

end 

if  mode3(k)<=0  &  mode3(k+l)>=0 
p=p+l; 

mode3_nodes(p)=pos(k) ; 
mode3_crossednode=mode3_crossednode+l ; 
elseif  mode3(k)>=0  &  mode3(k+l)<=0 
p=p+l ; 

mode3_nodes (p)=pos (k) ; 
mode3_crossednode=mode3_crossednode+l ; 

end 

end 


°/o 

subplot (4, 2, 4)  M4=V( : ,asize-6) ;  m4hz=abs (D (asize-6 , asize-6) ) /2/pi 
mode4=M4(l : 2 : 2*n_els) ;  raode4=real (mode4/max(mode4) ) ; 
plot (pos ,mode4) ; title (sprintf ( ’Mode4  703.2f  ’ ,m4hz)) 
set(gca,  ’XTick’  ,  [] ) 

°/0set  (gca,  ’Trick’  ,  [] ) 

i=0;  p=0;  mode4_crossedanti=0 ;  raode4_crossednode=0 ;  for  k=2:n_els-l 
if  mode4(k-l)<=mode4(k)  &  mode4(k+l)<=mode4(k) 
i=i+l ; 

mode4_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode4_crossedanti=mode4_crossedanti+l ; 

end 
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elseif  mode4(k-l)>=mode4(k)  &  mode4(k+l)>=mode4(k) 


i=i+l ; 

raode4_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode4_crossedanti=mode4_crossedanti+l ; 

end 

end 

if  mode4(k)<=0  &  mode4(k+l)>=0 
p=p+l; 

mode4_nodes (p) =pos (k) ; 
mode4_crossednode=mode4_crossednode+l ; 
elseif  mode4(k)>=0  &  mode4(k+l) <=0 
p=p+l ; 

raode4_nodes (p)=pos (k) ; 
mode4_crossednode=mode4_crossednode+l ; 

end 

end 


°/o 

subplot (4, 2, 5)  M5=V( : ,asize-8) ;  m5hz=abs (D (asize-8 , asize-8) ) /2/pi 
mode5=M5(l : 2 : 2*n_els) ;  mode5=real (mode5/max(mode5) ) ; 
plot  (pos  ,mode5) ;  title  (sprintf  ( ’Mode5  °/„3 .  If  ’  ,m5hz)) 
set(gca,  ’XTick’  ,  [] ) 

°/0set  (gca,  ’Trick’  ,  [] ) 

i=0;  p=0;  mode5_crossedanti=0 ;  raode5_crossednode=0 ;  for  k=2:n_els-l 
if  mode5(k-l)<=mode5(k)  &  mode5(k+l)<=mode5(k) 
i=i+l ; 
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mode5_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode5_crossedanti=mode5_crossedanti+l ; 

end 

elseif  mode5(k-l)>=mode5(k)  &  mode5(k+l)>=mode5(k) 
i=i+l ; 

raode5_maxcurve (i)=pos (k) ; 
if  k>=crackell  &  k<=crackel2 

mode5_crossedanti=mode5_crossedanti+l ; 

end 

end 

if  mode5(k)<=0  &  mode5(k+l)>=0 
p=p+l; 

mode5_nodes (p) =pos (k) ; 
mode5_crossednode=mode5_crossednode+l ; 
elseif  mode5 (k) >=0  &  mode5(k+l)<=0 
p=p+l ; 

raode5_nodes (p)=pos (k) ; 
mode5_crossednode=mode5_crossednode+l ; 

end 

end 


°/o 

subplot (4, 2, 6)  M6=V( : ,asize-10) ;  m6hz=abs(D(asize-10,asize-10))/2/pi 
mode6=M6 (1 : 2 : 2*n_els) ;  raode6=real (mode6/max(mode6) ) ; 
plot (pos ,mode6) ; title (sprintf ( ’Mode6  °/„3 .  If  ’  ,m6hz)) 
set(gca,  'XTick'  ,  [] ) 

°/0set  (gca,  ’YTick'  ,  [] ) 
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i=0;  p=0;  mode6_crossedanti=0 ;  raode6_crossednode=0 ;  for  k=2:n_els-l 
if  mode6(k-l)<=mode6(k)  &  mode6(k+l)<=mode6(k) 
i=i+l ; 

mode6_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

mode6_crossedanti=mode6_crossedanti+l ; 

end 

elseif  mode6(k-l)>=mode6(k)  &  mode6(k+l)>=mode6(k) 
i=i+l ; 

raode6_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode6_crossedanti=mode6_crossedanti+l ; 

end 

end 

if  mode6(k)<=0  &  mode6(k+l)>=0 
p=p+l; 

mode6_nodes (p) =pos (k) ; 
mode6_crossednode=mode6_crossednode+l ; 
elseif  mode6(k)>=0  &  mode6(k+l)<=0 
p=p+l ; 

mode6_nodes (p)=pos (k) ; 
mode6_crossednode=mode6_crossednode+l ; 

end 

end 


°/o 

subplot (4, 2, 7)  M7=V( : ,asize-12) ;  m7hz=abs(D(asize-12,asize-12))/2/pi 
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mode7=M7 (1:2: 2*n_els) ;  raode7=real (mode7/max(mode7) ) ; 
plot (pos ,mode7) ; title (sprintf ( ’Mode7  °/03.1f  ’,m7hz)) 
set(gca,  ’XTick’  ,  [] ) 

°/0set  (gca,  ’YTick'  ,  [] ) 

i=0;  p=0;  mode7_crossedanti=0 ;  raode7_crossednode=0 ;  for  k=2:n_els-l 
if  mode7(k-l)<=mode7(k)  &  mode7(k+l)<=mode7(k) 
i=i+l ; 

mode7_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode7_crossedanti=mode7_crossedanti+l ; 

end 

elseif  mode7(k-l)>=mode7(k)  &  mode7(k+l) >=mode7(k) 
i=i+l ; 

raode7_maxcurve (i)=pos (k) ; 
if  k>=crackell  &  k<=crackel2 

raode7_crossedanti=mode7_crossedanti+l ; 

end 

end 

if  mode7(k)<=0  &  mode7(k+l)>=0 
p=p+l; 

mode7_nodes (p) =pos (k) ; 
mode7_crossednode=mode7_crossednode+l ; 
elseif  mode7(k)>=0  &  mode7(k+l) <=0 
p=p+i; 

raode7_nodes (p)=pos (k) ; 
mode7_crossednode=mode7_crossednode+l ; 

end 
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end 


°/o 

subplot (4, 2, 8)  M8=V( : , asize-14) ;  m8hz=abs(D(asize-14,asize-14))/2/pi 
mode8=M8(l : 2 : 2*n_els) ;  raode8=real (mode8/max(mode8) ) ; 
plot (pos ,mode8) ; title (sprintf ( ’Mode8  °/„3 .  If  ’  ,m8hz)) 
set(gca, ’XTick’ , [] ) 

°/0set  (gca,  ’YTick'  ,  [] ) 

1=0 ;  p=0;  mode8_crossedanti=0 ;  mode8_crossednode=0 ;  for  k=2:n_els-l 
if  mode8(k-l)<=mode8(k)  &  mode8(k+l)<=mode8(k) 
i=i+l ; 

mode8_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

mode8_crossedanti=mode8_crossedanti+l ; 

end 

elseif  mode8(k-l)>=mode8(k)  &  mode8(k+l)>=mode8(k) 
i=i+l ; 

raode8_maxcurve(i)=pos(k) ; 
if  k>=crackell  &  k<=crackel2 

raode8_crossedanti=mode8_crossedanti+l ; 

end 

end 

if  mode8(k)<=0  &  mode8(k+l)>=0 
p=p+l; 

mode8_nodes (p) =pos (k) ; 
mode8_crossednode=mode8_crossednode+l ; 
elseif  mode8(k)>=0  &  mode8(k+l) <=0 
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p=p+l; 

mode8_nodes (p)=pos (k) ; 
mode8_crossednode=mode8_crossednode+l 


end 


end 


°/o 
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Appendix  C.  Beam  Drawings 


The  following  pages  display  diagrams  created  by  Aaron  Chrniel  of  the  eccen¬ 
trically  notched  beams.  The  dimensions  are  in  inches. 
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